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^ ABSTRACT 

CO Context. Searching for planets around stars with different masses probes the outcome of planetary formation for different initial con- 

ditions. The low-mass M dwarfs are also the most frequent stars in our Galaxy and potentially therefore, the most frequent planet 
hosts. 

JT"" Aims. This drives observations of a sample of 102 southern nearby M dwarfs, using a fraction of our guaranteed time on the 

Mh ESO/HARPS spectrograph. We observed 460 hours and gathered 1965 precise (~ 1 - 3 m/s) radial velocities, spanning the pe- 

PQ riod from Feb. 1 1th, 2003 to Apr. 1st 2009. 

Methods. This paper makes available the sample's time series, presents their precision and variability. We apply systematic searches 
for long-term trends, periodic signals and Keplerian orbits (from 1 to 4 planets). We analyze the subset of stars with detected signals 
and apply several diagnostics to discriminate whether the observed Doppler shifts are caused by stellar surface inhomogeneities or 
Q by the radial pull of orbiting planets. To prepare for the statistical view of our survey we also compute the limits on possible unseen 

?-H signals, and derive a first estimate of the frequency of planets orbiting M dwarfs. 

^ Results. We recover the planetary signals corresponding to 9 planets already announced by our group (Gl 176 b, Gl 581 b, c, d & e, 

^ G1674 b, G1433 b, G1667C b and G1667C c). We present radial velocities that confirm GJ 849 hosts a Jupiter-mass planet, plus 

' a long-term radial-velocity variation. We also present RVs that precise the planetary mass and period of Gl 832b. We detect long- 

term RV changes for G1367, G1680 and G1880 betraying yet unknown long-period companions. We identify candidate signals in 
vN the radial-velocity time series of 1 1 other M dwarfs. Spectral diagnostics and/or photometric observations demonstrate however that 

^ they are most probably caused by stellar surface inhomogeneities. Finally, we find our survey sensitive to few Earth-mass planets for 

periods up to several hundred days. We derive a first estimate of the occurrence of M-dwarf planets as a function of their minimum 
mass and orbital period. In particular, we find that giant planets (m sin / = 100 - 1, 000 M©) have a low frequency (e.g. / < 1% for 
O P = 1 - 10 d and / = 0.02^^ Jj^ for P = 10 - 100 d), whereas super-Earths (m sin / = 1 - 10 M©) are likely very abundant (/ = 0.36:j; 

^ for P = 1 - 10 d and / = 0.35+j; for P = 10 - 100 d). We also obtained t/© = 0.41+0;^^, the frequency of habitable planets orbiting 

M dwarfs (1 < m sin / < 10 M®). For the first time, 77© is a direct measure and not a number extrapolated from the statistic of more 
massive and/or shorter-period planets. 

I Key words. Stars: late-type - planetary systems, technique: radial- velocity 



k> 1. Introduction |Gatewood|1995| ), targeting M dwarfs has proven more success- 



ful since. 

J3 M dwarfs are the dominant stellar population of our Galaxy (e.g. 

IChabrier & Baraffe|2000| ). Compared to our Sun, they are cooler, At the forefront of planet discoveries, the radial-velocity 



smaller and lower-mass stars. These characteristics ease the de- (R^) technique was first to unveil a candidate giant planet or- 
tection of planets for many techniques and M dwarfs have there- biting an M dwarf. Thr ee years after the discovery of 5 IPeg b 
fore been included early in planet-search samples. While the first ( [Mayor & Queloz|1995|), the detection of a giant p lanet orbiting 



claimed detections (e.g. around Barnard's star - [van de Kamp| the M dwarf GJ 876 ( [Delfosse et al.|[T99 8i Marc y et al.||1998| ) 



[T963] ) were later found incorrect ( |Gatewood & Eichhom|| 19731 P^^^^^ ^^^^ ^ ^^^^^^ ^^^^^ P^^^^^^ ^^t^" 

ally one of a few tens of M stars monitored by radial-velocity 

Send offprint requests to: X. Bonfils surveys, and its detection made the early impression that giant 

e-mail: Xavier.Bonfils@obs.ujf-grenoble.fr planets could be common around late-type stars. Today, only 

* Based on observations made with the HARPS instrument on the 6 M dwarfs are known to host a planet with a minimum mass 

ESQ 3.6-m telescope at La Silla Observatory under programme ID > 0.5 Mj^p (see Table [T]) and the view has progressively shifted 

072.C-0488(E) toward a low rate of occurrence for giant planets, compared to 
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sun-like stars ([Bonfils et al.|2006t [Butler et al.|20Q4]|Endl et al 



2006| [Johnson et al.|20Q7[|Cum ming et al.|2Q08| ). 

Improvements of the RV technique has led to the discovery 
of lower-mass planets down to msin/ ^1.9 M® ( [Mayor et al.[ 
|2Q09 ). Below 25 M®, there are 8 known M-dwarf hosts and they 
altogether host 12 such low-mass exoplanets. Hence, despite an 
adverse detection bias, planets of low-mass already appear more 
frequent than giant planets ( [Bonfils et al.[[20Q7| ). Among them, 
GJ581d and GJ667Cc are noticeably interesting because they 
have msin/ < 10 M® and receive closely the amount of light 
received by Earth in our Solar System ( [Udry & Santos[[2QS7} 



Mayor et al. 2009, Delfosse et al., in prep.). Depending on their 



atmosphere (thickness, albedo and chemistry) liquid water may 
flow on thei r surface - the standard criterium to d efine a habit- 
able planet ( |Kasting et al.|1993[[Seisis et al.|200"7] ). 

The transit technique has also been successful in detect- 
ing two planets transiting an M dwarf. One is GJ 436 b, a 
Neptune-mass planet initially detected with Doppler measure- 
ments ( ,Butler et al. 2006t[Maness et al. 2007) and subsequently 
seen in transit ( [Gillon et al.[[2007b| ). Finding that GJ 436 b un- 
dergoes transits has enabled a wealth of detailed studies such 
as the determinations of the planet's true mass and radius and 
measurements of its efl'ective temperature and orbital ecce ntric- 
ity (Gillon et al. 2007a i Demory et al. 2007l[Deming et al.[2007) . 



Most recently, the Mearth pro ject, a search for transiting plan - 
ets dedicated to late M dwarfs ( [Nutzman & Charbonneau|2()08] ), 
has unveiled a ~ 6 M® planet transiting the nearby M4.5 dwarf 
GJ 1214 dCharbonneau et al.[ [2009^. Like GJ436b, it has a fa- 
vorable planetary to stellar radius ratio and is well suited to in- 
depth characterizations with current observatories. Both planets 
are considered Rosetta stones to the physics of low-mass planets. 

Anomalies in gravitational microlensing light curves can re- 
veal planetary systems kiloparsecs away from our Sun. Most fre- 
quently, the lenses are low-mass stars of masses < 0.6 M® and of 
spectral types M and K. Up to now, the technique has found 12 
planets in 1 1 planetary systems. Among those, 7 are giant planets 
and 5 fall in the domain of Neptunes and Super-Earths (Tab.[T]). 
The technique is mostly sensitive to planets at a few AUs from 
their host which, for M dwarfs, is far beyond the stellar habit- 
able zone. The microlensing technique probes a mass- separation 
domain complementary to the RV and transit techniques and has 
shown evidences that, at large separation s, low-mass planets out- 
number giant planets ( [Gould et al.|2006| ). 

Ground-based astrometry applied to planet searches has been 
cursed by false positives, of which Van de Kamp's attempts 
around Barnards star are probably the most fam ous examples 
(van de Kamp 1963[ Gatewood & Eichhom[[1973( ). Fifty years 
ago, Van de Kamp first claimed that a 1.6 Mj^p planet orbits 
Barnard's star every 24 years. Over the following decades, he 
continued to argue for a planetary system around the star (van de| 
[Kamp|p^82]), despite growing evidence of systematic s in the 
data (e.g. [Gatewood & Eichhorn|1973| |Hershey'1973). Radial- 
velocity and astrometric data have now completely excluded 
the van de Kamp pla nets (Gat ewood|1995[[Kurster et al.||2003[ 
[Benedict et al.[[2002| ), but Barnard's star has been far from the 
only target with false astrometric detections. 

Nevertheless, astrometry has proven useful to confirmed the 
planetary nature of a few radial- velocity detections, and to re- 
move from the planet sample the low-mass stars s een with an 
unfavorable inclination (e.g. Halbwachs et al.|2()00| ). Moreover, 
thanks to HST/FGS astrometric observations, GJ 876b has been 



the second exoplanet with a true mass determination ( [Benedict 
et al.|2002) ), soon after the detection of the transiting planet HD 



209458 b ([Charbonneau et"aL][2Q00t [Mazeh et al.|[2000l [Henry 
let al.|2000 V 

To complete the view of planetary-mass objects formed at 
the lower end of the main sequence, let us mention the ~ 5 Mjup 
companion detected with RV measurements around the young 
brown dwarf Cha Ha 8 ( [Joergens & Muller[ |2007| ) and the 
~ 5 Mjup companion imaged around another young brown dwarf 
2M1207 ( Cha uvin et a l. 2004). The protoplanetary disks of both 
brown dwarfs were most likely not massive enough to form 
such massive objects, since observations show that protoplan- 
etary disk masses scale at most linearly with the mass of the star. 
Both 2M1207b and Cha Ha 8b therefore probably formed like 
stars rather than in protoplanetary disks. 

Table [T] lists the known M-dwarf hosts and their procession 
of planets. For each planet, it gives the basic characteristics and 
a reference to the discovery papers. In total, there are 35 planets 
in 28 planetary systems. 

Planets orbiting M dwarfs formed in a diff'erent environment 
than those around solar-type stars, and therefore reflect a difl'er- 
ent outcome of the planetary formation mechanism. The mass 
of the proto-planetary disk, its temperature and density profiles, 
gravity, the gas-dissipation timescale, etc... all change with stel- 
lar mass (e.g. Ida & Lin,200 5). For the construction of the Harps 
spectrograph for ESO, our consortium has been granted 500 ob- 
serving nights with the instrument spread over 6 years. We chose 
to dedicate 10% of that guaranteed time to characterize the plan- 
etary population for stars with masses < 0.6 M®. 

This paper reports on our 6-year RV search for planets 
around M dwarfs and the outline is as follow. We first de- 
scribe our sample (Sect. [2]) and present the RV dataset collected 
(Sect. [3]). We next perform a systematic analysis for variability, 
long-term trend and periodic signals (Sect. [4]). We close that sec- 
tion by an automated Keplerian multi-planet search. For all sig- 
nals detected with enough confidence, we apply a suite of diag- 
nostics to disentangle Doppler shifts caused by bona fide planets 
from Doppler shifts caused by stellar surface inhomogeneities 
(Sect. [5]). Sect. 6 presents the detection limits for individual stars 
of our sample and Sect. 7 pools them together to estimate both 
the survey sensitivity and the frequency of planets orbiting M 
dwarfs. Sect. 8 summarizes our results and presents our conclu- 
sions. 



2. Sample 

Our search for planets orbiting M dwarfs originates from RV ob- 
servations started in 1995 with the 1.93m/ELODiE spectrograph 
(CNRS/OHP, France). This former program aimed at determin- 
ing the st ellar multiplicity for very-low-mass stars (Delf osse] 
et al.|1999) as well as the detection of their most massive plan- 
ets ( [Delfosse et al.1[1998| ). In 1998, we started to run a similar 
program in the southern hemisphere, with the 1.52m/FEROS spec- 
trograph (La Silla, Chile). With Feros's early settings, we were 
unsuccessful to improve on its nominal precision of ~ 50 m/s. 
Nevertheless, we benefitted from these observations to start with 
a better sample on Harps, cleaned from spectroscopic binaries 
and fast rotators for which precision radial- velocity is more dif- 
ficult. 

Our Harps sample corresponds to a volume limited of M 
dwarfs closer than 11 pc, with a declination d < -h20°, brighter 
than V = 14 mag and with a projected rotational velocity 
vsin/ < 6.5 km/s. We removed known spectroscopic binaries 
and visual pairs with separation < 5'' (to avoid light contam- 
ination from the unwanted component). We have however en- 
countered a few spectroscopic binaries and fast rotators which 
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Table 1. Known exoplanets orbiting M dwarfs and their basic parameters 



Radial-Velocity detections 



Star 


PL 




sin 


i 


P 


a 


e 


ref. 


ref. 


in 






[Me] 


[Mjup] 


n 


[d] 


[AU] 




discovery* 


param.* 


sample 


GJ176 


b 


8.4 


0.026 


- 


8.7(8) 


0.07 


0. (fixed) 


25 


25 


y 


GJ317 


b 


380 


1.2 


- 


69(2) 


0.95 


0.19(006) 


17 


17 


n 


GJ433 


b 


6.41 


0.0202 


- 


7.36(5) 


0.06 


0. (fixed) 


37 


37 


y 


GJ581 


b 


15.7 


0.0492 


>40. 


5.368(7) 


0.04 


0. (fixed) 


8 


23 


y 




c 


5.4 


0.017 


>40. 


12.9(3) 


0.07 


0.17 (0.07) 


13 


23 


y 




d 


7.1 


0.022 


>40. 


66.(8) 


0.22 


0.38 (0.09) 


13,23 


23 


y 




e 


1.9 


0.0060 


>40. 


3.1(5) 


0.03 


0. (fixed) 


23 


23 


y 


GJ649 


b 


106 


0.333 


- 


59(8) 


1.14 


0.30 (0.08) 


28 


28 


n 


GJ667C 


b 


6.0 


0.019 


- 


7.20(3) 


0.05 


0. (fixed) 


37 


37 


y 




c 


3.9 


0.012 


- 


28.1(5) 


0.28 


0. (fixed) 


37 


37 


y 


GJ674 


b 


11 


0.034 


- 


4.6(9) 


0.04 


0.20 (0.02) 


12 


12 


y 


GJ676A 


b 


1300 


4.0 


- 


98(9) 


1.61 


0.29 (0.01) 


36 


36 


n 


GJ832 


b 


200 


0.64 




3(416) 


3.4 


0.12(0.11) 


20 


39 


y 


GJ849 


b 


310 


0.99 




18(52) 


2.35 


0.04 (0.02) 


11 


39 


y 


GJ876 


b 


839 


2.64 


48.(9) 


61.0(7) 


0.211 


0.029 (0.001) 


1,2 


29 


y 




c 


180 


0.83 


4(8) 


30.2(6) 


0.132 


0.266 (0.003) 


3 


29 


y 




d 


6.3 


0.020 


50 (fix.) 


1.9378(5) 


0.021 


0.139 (0.032) 


6 


29 


y 


GJ3634 


b 


6.6 


0.021 




2.645(6) 


0.028 


0.08 (0.09) 


36 


36 


n 


HIP 12961 


b 


110 


0.35 




57.4(3) 


0.13 


0.16(0.03) 


35 


35 


n 


HIP 57050 


b 


40 


0.3 




41.(4) 


0.16 


0.3 (0.1) 


34 


34 


n 


HIP 79431 


b 


350 


1.1 




111.(7) 


0.36 


0.29 (0.02) 


30 


30 


n 



Transit detections 



Star 


PI. 




ml 

[Mjup] 




P 


a 


e 


Rp 


ref. 


ref. 


in 






[Me] 


n 


[d] 


[AU] 






discovery* 


param.* 


sample 


GJ436 


b 


22.6 


0.0711 


85.(9) 


2.643(9) 


0.029 


0.14(0.01) 


4.(2) 


5,14 


15, 16 


n 


GJ1214 


b 


6.5 


0.020 


88.(6) 


1.5803(9) 


0.014(3) 


<0.27 


2.(7) 


26 


26 


n 



Microlensing detections 



Star 


PI. 






a 




ref 


ref. 






[Me] 


[Mjup] 


[AU] 


[Mo] 


discovery* 


param.* 


OGLE235-MOA53 


b 


830-- 




4.3+2-^ 


0.67+0.14 


4 


22 


MOA-2007-BLG-192-L 


b 






^•'^'^-0.06 


0.084!^:0}^ 


19 


32 


MOA-2007-BLG-400-L 
OGLE-2007-BLG-368-L 


b 
b 


260:i60 


0.83!g 
0.23^0.05 


72+0-3^ 1 6 5+^-2 ** 

o + 1.4 
•^••^-0.8 

1.25+0.10 


0.30!^j^ 
0.67+0.14 


24 
31 


24 
31 


MOA-2008-BLG-310-L 


b 


74+17 


27 


27 


OGLE-06-109L 


b 


226+25 


0.711+0.079 


2.3+0.2 


0.50+0.05 


18 


22 




c 


86+10 


0.27+0.03 


4.6+0.5 




18 


22 


OGLE-05-169L 


b 


13:^ 


QQ4^+0.013 


3 2+^-^ 


Q 49+0.14 

22+^-^^ 

^•^^-0.11 


10 


22 


OGLE-05-390L 


b 


5.5+^-^ 


017+^-^'^ 

^•^^ ^-0.008 


26;1:? 


9 


124 


OGLE-05-071L 


b 


1200+100 


3.8+0.4 


2.1+0.1 13.6+0.2 ** 


0.46+0.04 


7 


21 


MOA-2009-BLG-319-L 


b 


50!^^ 


0.2+0.1 


24+1.2 


0.38:^-34 


33 


33 


MOA-2009-BLG-387-L 


b 


830^300 




1 o+o:^ 

^•^-0.7 


38 


38*** 



* The true mass (m^) is reported for GJ 876 b, c, for the transiting planets GJ436b and GJ 1214b and for all microlensing detections. The masses 
given for G/876J assumes a 50° orbital inclination. We give minimum masses GJ581b, c, d and e, and dynamical consideration restrict coplanar 
systems to / > 40°. Usually, uncertainties in planetary masses do not include the stellar mass uncertainty. 

** degenerated solution 

instead of Icr uncertainties, we quote 90 % confidence intervals fromlBatista et al. |(|201l|) 



;(5) 



Butler et al. 



(120041; (6) 



Rivera et al 



(|2005 



* (1) Delfosse et al. ( 1998 ); (2) Ma rcy et al. |(1998 ); (3) Marcy et al. (2001 ); (4) Bond eTaL||2004l . 

(7)^Udalski et al. (2005 ); (8) Bonfils et al. (2005 ); (9) Beaulieu et al. (2006); (10) Gould et al. (2006); (ll)"Butler et al."b00 6|); (12)|Bo nfils et al.| 
([2007); (13) Udry & Santos (2007); (14) Gillon et al. (2007b); (15) Gillon et al. (2007a); (16) Demory et al. (2007); (17) Johnson et al. (2007); 
(18) Gaudi et al. (2008 ); (19) Bennett et al. (2008 ); (20) Bailey et al. (2009 ); (21) Dong et al. (2009b); (22) Bennett (2009 ); (23) Mayor et al] 
([2009 ); (24) Dong et al. (2009a); (25) Forveille et al. ( 2009 ); (26) Charbonneau et al. (2009 ); (27) Janczak et al. ( 2010); (28) Johnson et al. ( 20T0k 
(29) Correia et al. (2010); (30) Apps et al. (2010); (31)|Sumi et al.| ( |2010| ); (32)|Kubas et al. (2010); (33) Miyake et al. (2010); (34) Haghigh ipour] 
|etal.. ( ,2010j (3 5) . Forveille et al., ( .2011j ; (36) Bonfils et al.| ( |2011| ); (37) Delfosse (2011, in prep.); (38) ,Batista et al.,(,2011j ; (39) This paper. 



were not known before our observations. We list them in Table [2] 
and discard them from the sample presented here. Note that 
we also dismiss GJ 1001, Gl 452.1 and LHS 3836. The first two 
stars were initially counted in our volume limited sample (with 
TT = 103 and 96 mas, respectively - jOliese & JahreiB|1991| ) and 



had their parallax revised since, now placing them beyond 1 1 pc 
(with 7T = 76.86 + 3.97 and 88.3 ± 3.7 mas - [Henry et al.|2006| 
Smart et al.|2010|). LHS 3836 was initially included based on its 



V magnitude in Gliese & JahreiB[ ( [T991 )'s catalog but our first 
measurements were indicative of a much lower brightness. 
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Table [3] lists the 102 stars selected for the sample. Their co- 
ordinates, proper motions and parallaxes are primarily retrieved 
from the revised Hipparcos catalog ( van Leeuwen 2 007 ). A frac- 
tion of the parallaxes, unavai lable in the Hipparcos database , 

( p^95] ), [Henry eTaLl ( [20Q6| ), 



were obtained from|^van Altena et al. 



Reid et al.| ([T995J and the 4th Cata 



og of Nearby Stars (CNS4 

- Jahreiss, priv. comm.). V-band magnitudes are taken from 
Simbad and infrared J- and K-band magnitudes from 2MASS 



( |Cutrietal.|20Q3 ). We used the empirical |Delfosse et al.| ( |2000l )'s 
mass-luminosity relationship together with parallaxes and K- 
band photometry to compute the mass of each star. Infrared K- 
band photometry a nd (J- K) colors are converted to luminosities 
with [Leggett et al. ( 2001| )'s bolometric correction. 

We also indicate in Table [3] the inner and outer limits for the 
distance of the Habitable Zone using the recent- Venus and early- 
Mars criterions, respectively, and Eq. (2) and (3) from Selsil| 
|et al.| ( [2Q07| ). The boundaries of the Habitable Zone are uncertain 
and depend on the planet's atmospheric composition. Extra-solar 
planets found close to these edges have therefore to meet more 
stringent conditions to be inhabitable. For more detailed consid- 
erations, we refer the reader to more comprehensive models (e.g. 
ISelsis etal.|2QQ7l ). 

Our sample is composed of the closest neighbors to the Sun. 
Nearby stars tend to have large proper motions and the projection 
of their velocity vector may change over time, up to few m/s/yr 
( [Schlesinger] 1 9 1 7 [ Kurster et al.^2003 ) . We therefore report the 
value of their secular acceleration in Table O 

To portray our sample, we show its V-mag. and mass dis- 
tributions in Fig. [T] For both distributions, the average (resp. 
median) value is plotted with a vertical straight (resp. dashed) 
line. The magnitudes and masses of planet hosts are also marked 
with vertical ticks on top of the histograms. The target bright- 
ness spans V=7.3 to 14 mag. with a mean (resp. median) value 
of 11.25 mag. (resp. 11.43 mag.). The stellar mass ranges from 
0.09 to 0.60 Mq with an average (resp. median) value of 0.30 
Mo (resp. 0.27 M©). The smaller count seen in the 0.35-0.40 M© 
bin is unexplained but from statistical fluctuations. Interestingly, 
one can note that our sample covers a factor of ~ 6 in stellar mass, 
while the mass step between our typical M dwarf (~ 0.27 M©) 
and the typical Sun-like star (~ 1 M©) corresponds to a factor 
of less than 4. This means that planetary formation processes 
depending on stellar mass could lead to larger observable difl'er- 
ences across our sample than between our M-dwarf sample and 
Sun-like stars. 

There are overlaps between our sample and others that sim- 
ilarly targets M dwarfs to search for planets. Among them, we 
found published RV time-series for Gl 1, Gl 17 6, G1229 , G1357, 
G1551, G1682, G1699, G1846 and G1849 in |Endl et al.| ( |2Q06l 
hereafter E06), for Gl 1, G1 229, G1357, G1433, G1551, G1682, 
Gl 699, Gl 846 and Gl 849 in |Zechmeister et al.| (2009 , hereafter 
Z09), and few others in detection papers, as for Gl 176 (|Endl 



et al. 



2008)), Gl 832 (B ailey et al.|20 09 ) and Gl 849 (B utler et al. 



2006| ). When possible, we compare our results to these time se- 
ries and, for completeness, additional comparison is given in 
Appendix [A| 



3. Observations 

To gather RV observations for the sample described above we 
used the Harps instrument ( Mayor et al.||2003t |Pepe et aL] 
2004]), a spectrograph fiber fed with the ESO/3.6-m telescope 




E 15 - 



V [mag.] 




0.1 0.2 



0.3 0.4 
Mass [M^] 



(La Silla, Chile). It covers the 3800-6800 A wavelength do- 
main, has a resolution R ~ 115 000 and an overall throughput 



Fig. 1. Sample distributions for V magnitudes and stellar 
masses. The vertical dashed and plain lines locate the median 
and averaged values, respectively. The small ticks are explained 
in Sect.E 



(instrument-h telescope) greater than 6%. It is enclosed in a vac- 
uum vessel, is pressure controlled to ±0.01 mbar and thermally 
controlled to ±0.01 K. That ensures a minimum instrumental 
shift on the position of the spectrum image on the CCD or, prac- 
tically, a RV drift < 0.5 m/s/night. To reference the origin of 
the instrumental RV drift, the fiber can be illuminated with a 
Thorium-Argon (ThAr) lamp at any time. Harps also off'ers a 
second fiber that can be illuminated with ThAr light simulta- 
neously while the scientific fiber receives star light. This mode 
avoids the need to record frequent calibrations between scien- 
tific exposures and can correct the small instrumental drift that 
occurs while the stellar spectrum is recorded. Since the instru- 
mental drift during the night is small, this mode is only used 
when a sub-m/s precision is required. Our observational strat- 
egy for M dwarfs aims to achieve a precision of ~ 1 m/s per 
exposure for the brightest targets. We chose therefore not to use 
the second fiber and relied on a single calibration done before the 
beginning of the night. As the science and calibration spectral or- 
ders are interlaced, avoiding the second fiber eludes light cross- 
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5 



Name 


a (2000) 


6 (2000) 


V [mag] 


Comment 


L 225-57 


02:34:21 


-53:05:35 


7.3 


SB2 


CD-44-836B 


02:45:14 


-43:44:06 


12.7 


Fast Rotator /SB2 


LHS 1610 


03:52:42 


-hl7:01:06 


13.7 


SB2 


LHS6167 


09:15:36 


-10:35:47 


14.7 


SB2^ 


G 161-71 


09:44:55 


-12:20:53 


13.8 


Fast Rotator 


GJ1154A 


12:14:17 


-h00:37:25 


13.7 


Unresolved SB 2* 


LHS 3056 


15:19:12 


-12:45:06 


12.8 


SB 


L 43-72 


18:11:15 


-78:59:17 


12.5 


Unresolved SB2* 


LTT7434 


18:45:57 


-28:55:53 


12.6 


Unresolved SB2* 


G1867B 


22:38:45 


-20:36:47 


11.4 


Fast Rotator or unresolved SB 2 



" previously detected by |Montagnier et al.| ( |2006 1 
X SB with variable spectral-line width 



contamination between science and calibration spectra. This can 
be a source of noise for the blue-most spectral orders, where we 
only reach low signal-to-noise for M dwarfs. In particular, we 
were interested in clean Ca ii H&K lines because they are a use- 
ful diagnostic of stellar activity (see Sect. [5]). 

From the first light on Harps on February, 1 1th 2003 to the 
end of our guaranteed time program on April, 1st 2009, we have 
recorded 1965 spectra for the M-dwarf sample, for a total of 
460 h of integration time. 

We computed RV by cross-correlating the spectra with a nu- 
merical weig hted mask following [Baranne et al.| ( |1996] ) and |Pepe 
et al.| ( |20()2| ). Our numerical mask was generated from the addi- 
tion of several exposures taken on Gl 877, a bright M2 star of 
our sample. Co-addition of spectra requires knowing their rela- 
tive Doppler shifts. We computed RVs for Gl 877 spectra with 
a first iteration of the template and re-built the template more 
precisely. We obtained convergence just with a few iterations. 
The numerical mask counts almost 10,000 lines and most of the 
Doppler information in the spectrum. No binning is done. 

The RV uncertainties were evaluated from the Doppler con- 
tent of individual spectra, using the linear approximation of a 
Doppler shift (Eq. 12 - Bouchy et al. 2001 ). This formula gives 
more weight to spectral elements with higher derivative because 
they are more sensitive to phase shifts and contribute more to 
the total Doppler content. It is important to note that we do not 
sum the Doppler content of individual spectral elements over 
the whole spectrum. The derivative of the spectrum has a higher 
variability against noise than the spectrum itself and, for low 
signal-to-noise ratio, doing so would over-estimate the RV pre- 
cision. Instead, to mitigate the eff'ect and compute more realis- 
tic uncertainties, we applied the formula directly to the cross- 
correlation profile, which has ~ 30 times higher signal-to-noise 
than t he individual spectral lines (see appendix A in ,Boisse et al.| 
2Q1Q| ). To account for the imperfect guiding (-30 cm/s) and 
wavelength calibration (~50 cm/s) we quadratically added 60 
cm/s to the Doppler uncertainty. 

As a trade-off' between exposure time and precision we chose 
to fix the integration time to 900 s for all observations. We ob- 
tained a precisiorj^cr^ ~ 80 cm/s from V^^^ = 7-10 stars and 
(Ti ~ 2.5^^^-^^/^ m/s for V^^^ = 10 - 14. Our internal errors cr, 
(composed of photon noise + instrumental errors) are shown in 
Fig. |2] where we report, for all stars with more than 6 measure- 
ments, the mean ct/ (blue filled circle) as a function of the star's 
magnitude, with error bars corresponding to ct/'s dispersion. For 



^ As opposed to our precision, our measurement accuracy is poor. 
Absolute radial velocities given in this paper may not be accurate to + 
1 km/s. 




Fig. 2. Internal (cr/) and external (cTg) errors as a function of 
V-band magnitudes, for stars with 6 or more measurements. 



comparison. Fig. [2] shows observed dispersions CTe for all stars 
with more than 6 measurements (black squares, changed to tri- 
angles for clipped values). The cr^ are the observed weighted 
r.m.s. and are related to the;^^ value : 



, Y,{Ry-<Ry>Ylo-' 
(t: = — 



Zi/^ 



(1) 



We discuss the diff'erence between internal and external errors in 
Sectiongl 

Our RV time series are given in the Solar System barycentric 
reference frame and can be retrieved online from CDS. Prior to 
their analysis, the time series are corrected from the secular RV 
changes reported in Tab. [3] We also show the time series in Fig.|3] 
(only available on-line), after subtraction of half the min-hmax 
value for a better readability. 



4. Data analysis 

Several planet- search programs have presented a statistical anal- 



ysis for their survey (e .g. Murdoch et al. 1993; Walker e t al 
ICumming qV^\999 , 2008 ; Endl et al. 2002 ; Zec hmeister 



1995 



et al. 



2Q09| ). Often, statistical tests are applied to the time se- 



ries in order to appraise the significance of trends or variability. 
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Then, the time series are searched for periodicities and, if a sig- 
nificant periodicity is found, the corresponding period is used as 
a starting point for a Keplerian fit. Again, statistical tests are ap- 
plied to decide whether a sinusoidal or a Keplerian model is a 
good description for the time series. In this section, we follow 
the same strategy and add a heuristic method based on genetic 
algorithms to the systematic search for Keplerian signals. 



4.1. Excess variability 

Once we have computed the RVs, the Doppler uncertainties (ct/) 
and the RV dispersion (cr^) for each star, the first step is to test 
our time series for variability. In Fig.|2j we have reported both ct/ 
and (Tg as a function of stellar magnitudes for all stars with more 
than 6 measurements. Apart from one star (Gl 447, which has 6 
RVs only), we have cr^ > cri. The cr^ have a lower envelope that 
matches the ct/'s envelope for V"^^^ = 10-14 and that is slightly 
above (~ 2 m/s) in the brighter range y^^^ = 7-10. 

To test whether the observed RVs vary in excess of our inter- 
nal errors, we first compare the expected variance < cr^ > (the 
mean internal error) to the measured variance cr^ by applying a 
F-test, which gives a probability P(F) to the F-value F = 



data sets with xlirtuai ^ ^liope gi^^^ the false-alarm probabil- 
ity (FAP) for the slope model. For all stars, but the ones with 6 
measurements or less, we generate 1,000 virtual time series. And 
because that method probes FAPs greater than 0(1 /N\) only, we 
limit the number of trials to Nl for stars with fewer measure- 
ments. 

Table |4] gives the slope coefficient a as well as P(F) and 
FAP values for all time series. The P(F) and the FAP values 
are reported in boldface when below a threshold of 1%, to point 
high confidence level in the slope model. 

Among our sample and according to the FAP values, we 
find that 15 stars have time series better described by a slope 
than a simple constant. They are Gil, LP771-95A, G1205, 
G1341, G1382, G1413.1, G1618A, G1667C, G1680, G1699, 
G1701, G1752A, G1832, G1849 and G1880. We also see that, 
while LP 77 1-95 A, Gl 367, G1618A, Gl 680 and Gl 880 displays 
smooth RV drifts, the other 10 stars seems to obey a more com- 
plex variability. According to P(F) values, we find that the same 
15 stars plus 8 more have a significant chi squared improvements 
when we fit a slope. They are G154.1, G1250B, G1273, G1367, 
G1433, G1551, G1674 and G1887. 



(e.g. |Zechmeister et al.|2009| ). As another test for variability, we 4.3. Periodicity 



also compute the xionstant^^^ the constant model and Pixlonstam)^ 
the probability of hsLwing xlonstam given the cr^. For both the F- 
test and the;^^-test, a low probability means that photon noise, 
calibration uncertainty and guiding errors are not sufficient to ex- 
plain the observed variability. In such a case, one has to invoke 
an additional source for variability, from unaccounted noise to 
planetary companions. 

We report ct/, o-g, P(F), and P(x^) in Table |4] and change 
the P{F) and the Pix^) values to boldface when smaller than 
1%, i.e. when they indicate a confidence level for variability 
higher than 99%. Using this criterion, the F probabilities (resp. 
the;^^ probabilities) indicate that 45% (resp. 63%) of our sam- 
ple displays an excess variability. When focusing on stars with 
ymag = 7 _ IQ, all stars but two are found more variable than 
expected according to P{F), and all stars according to Pix^). 
The reason is that our external error never reaches the ~70 cm/s 
threshold estimated for the brighter range of our sample (but 
rather 1.5-2 m/s), and is dominated by photon noise for a third 
of the sample in the y^^^ = 10-14 range. 



4.2. Trends 

Next, we examine the time series for possible trends. They may 
correspond to incomplete orbits and betray the presence of long- 
period companions. For each star, we adjust a slope a (RV = 
at+P) to the RV data and evaluate the x^siope ^^^^^ that model. 

To know whether the slope is a better description to the data 
than the no slope model, we use two statistical tests. First, we use 



the F-test to gauge whether a lower x^i^p^ 



compared to;^^ 



a suflftcient improvement to justify an additional free parameter 
(2 for the slope model against 1 for the constant model). 

Also, because the F-test statistics is ill-behaved for non nor- 
mally distributed uncertainties, we use a less-sensitive test based 
on bootstrap randomization. That time, we generate virtual RV 
time series by shuffling the original data, i.e., we keep the same 
observing dates and attribute to each date a measurement ran- 
domly chosen among the other dates, without repeating twice 
a given measurement. On each virtual time series, we adjust a 
slope and compute its xlirtuai ^^1^^- Then, the fraction of virtual 



In our search for planets, variability selects the stars to focus 
on and trends reveal yet uncompleted orbits. Our next step in 
the search for planetary candidates is to look for periodic sig- 
nals. The classical diagnostic for periodic coherent signals in 
unevenly spaced time series is the Lomb-Scargle periodogram 
(Lomb 1976; Scargle 1982 ) or, to account for the unknown sys- 
tem's mean velocity, its generalized version, the ffoating-mean 
periodogram ( Gumming et al.^19 99 ). 

We therefore compute generalized Lomb-Scargle peri- 
odograms for all our time series with at least 6 measurem ents. 
We follow the formalism developed in Zech meister & Kurster| 
( 2009 ), and choose a power normalization where 1.0 means that 
a sinusoidal fit is perfect (x^ = 0) and where 0.0 means that a 
sinusoidal fit does not improve over a constant model. We cal- 
culate false-alarm probabilities very similarly to our trend anal- 
ysis, with a bootstrap randomization (preferred to a F-test). We 
create 1,000 virtual time series by shuffling the original data set. 
For each individual data set we compute a periodogram and lo- 
cate the highest power. Considering all periodograms from all 
virtual data sets, we compute the distribution of the power max- 
ima. The power value that is above 99% of all trial powers is then 
considered as the 1 % FAP. More generally, we attribute a FAP 
to the maximum power value found in the original data set by 
counting the fraction of the simulated power maxima that have a 
larger value. 

We show all periodograms in Fig. [16] (only available in elec- 
tronic format). The periods corresponding to the highest power, 
the corresponding signal semi-amplitude, the;^^, as well as the 
associated FAPs are reported in Table [5] As previously, we 
boldface the significant signal detections, i.e. the periods with 
a power excess that has a FAP < 1%. 

In ideal cases, long-term variability only aff'ects the long- 
period range of the periodogram. However, the sparse sampling 
of our time series can cause power excess from a long-period 
signal to leak to shorter periods. Removing the long-term signif- 
icant trends cleans the periodogram and may reveal periodic sig- 
nal at shorter periods. Thus, we also compute the periodogram 
for the time series after removing the slope adjusted in the pre- 
vious section, as well as corresponding false-alarm probabilities 
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(FAP2 in Table [5]). We record a noticeable change for G1618A 
and Gl 680 only. No power excess remains in their periodograms 
after subtraction of the slope, meaning their RV variation are 
mostly linear drifts. 

Our periodicity analysis identifies 14 stars with power excess 
with FAP2 < 1% are : Gl 176, G205, G1273, G1358, G1388, 
G1433, G1479, G1581, G1667C, G1674, G1832, Gl 846, Gl 849 
and Gl 876. The star Gl 887 additionally has FAP approaching 
our 1% threshold. 

After identifying possible periodicities, we use the candi- 
date periods as starting values for Keplerian fits. We search the 
residuals for periodic signals again, computing periodograms 
(see Fig. 17 in the online material) and locating power excesses 
with FAP < 1%. We find a probable second periodic signal for 
Gl 176, G1205, G1581, G1674 and G1876 (Table[6]). We notice 
that Gl 667C, Gl 832 and Gl 846 have FAP approaching our 1% 
threshold. Although not shown here, a third iteration only reveals 
coherent signals for Gl 581 and, possibly, G1667C. A fourth iter- 
ation found signal for Gl 581 only. 

4.4. Keplerian analysis 

Even in its generalized form, the Lomb-Scargle periodogram is 
optimal for sinusoidal signals only. Eccentric orbits spread the 
spectral power on harmonics, decreasing the power measured at 
the fundamental period, and fitting a Keplerian function is an ob- 
vious improvement. Compared to a periodogram search, it can 
detect planets with high eccentricity against a higher noise. It 
is often not used because it is non-linear in some of the orbital 
parameters. Traditional non-linear minimizations can only con- 
verge to a solution close to a starting guess, possibly far from the 
global optimum. Applying non-linear minimization from many 
starting guesses becomes quickly impracticable when the num- 
ber of planets increases. Keplerian functions are indeed transcen- 
dent and evaluating the radial velocity at a given time therefore 
need an integration which is computationally expensive. Finally, 
the sequential approach outlined above requires a higher signal- 
to-noise for systems with several planets on commensurable or- 
bits (i.e. with RV pulls of similar amplitudes). 

To work around these shortcomings, we make use of a hybrid 
method based on both on a fast non-linear algorithm (Levenberg- 
Marquardt) and genetic operators (breeding, mutations, cross- 
over). The algorithm has been integrated by one of us (D.S.) 
in an orbit analysis software named Yorbit. We give a brief 
overview of this software here, but defer its detail description 
to a future paper (Segransan et al. in prep.). While a popula- 
tion of typically 1, 000 solutions evolves, the top layer of Yorbit 
evaluates the performances of the diff'erent minimization meth- 
ods and tunes their time allocation in real time. Genetic algo- 
rithms efficiently explore the parameter space on large scales. 
They thus score well and are given more CPU time during the 
early phase of the minimization. Once promising solutions are 
found, non-linear methods converge much more efficiently to- 
ward the local optima. Hence, when new solutions arise outside 
of local minima, the non-linear methods are given more CPU 
time. This heuristic approach has proved very efficient and the 
solution to multi-planet systems is found in few minutes only, 
on common desktop computers. 

We search for planets using Yorbit on all stars with more 
than 12 measurements, neglecting planet-planet interactions at 
this point. We scaled the complexity of the tested models with 
the number of measurements. Although in principle 5N-hl RVs 
are enough to obtain a Keplerian fit to an N-planet system, we 
wish to minimize the number of spurious solutions and arbitrary 



require at least 12 RVs per planet in the model. Hence we use a 
1 -planet model for 12 RVs, both a 1 -planet and a 2-planet model 
for stars with more than 24 RVs, and 4 diff'erent models (1-4 
planets) for stars with more than 48 RVs. To complement those 
models we also use the same suite of models with the addition of 
a linear drift. We allow Yorbit to run for 150 seconds per planet 
in the model. 

Just as for the periodicity analysis, evaluating the credence of 
the model is essential. The of solutions necessarily improves 
for more complex models as the number of degree of freedom 
increases, and we want to evaluate whether this improvement is 
statistically significant or occurs by chance. As previously, we 
generate virtual datasets by shuffling the original data and re- 
taining the dates. We create 1,000 virtual datasets and refit all 
tested models 1 ,000 times with Yorbit. For each star and model, 
we obtain 1,000;^^ values and count how many are below the;^^ 
measured on the original data. This gives the FAP for that model, 
compared to the no-planet hypothesis. A model is considered to 
improve significantly when less than 1% of the virtual trials 
give as low a;^^. 

Once we find a significant model, more complex models are 
then evaluated against that model, and not against simpler mod- 
els anymore. We consider signals in that model are detected (i.e., 
for instance, we assume the system is composed of 2 planets if 
that model is a 2-planet model). To generate the virtual datasets 
we use the residuals around the best solution for the new ref- 
erence model (i.e., in our example, the residuals around the 2- 
planet model). Shuffling the residuals (and retaining the dates), 
we again create 1,000 virtual datasets and fit the more complex 
models using Yorbit. How many from these virtual trials are 
lower than the best;^^ obtained on the actual RVs then gives the 
FAP for the complex models, compared to the simpler model. 

We report the parameters of the best solutions for each star 
and model in TablejT] In Tables [8]and[9j we report the FAP of the 
models, compared to selected simpler models. In both tables, we 
boldface the models with FAP<1%. Of the 43 stars with more 
than 12 measurements, 19 have RVs well modeled by 1 or more 
planets. Among them, we recover all stars with probable RV pe- 
riodicity except G1680 that has less than 12 measurements and 
was not tested. In the following section, we discuss the interpre- 
tation of these candidate signals. 



5. Interpretation 

The above analysis reveal Keplerian signals in our time series, 
but Doppler shifts do not always correspond to planets. To vet 
a RV variability against stellar activity, we make use of several 
diagnostics. The shape of the cross-correlation function (CCF) 
in particular is often informative. While its bary center measures 
the radial- velocity, its bisector inverse slope (BIS 



Queloz et al. 



2001 ) and its full width half maximum (FWHM) can diagnose 
inhomogeneities of the stellar surface. Alternatively, spectral in- 
dices built on Ca 11 H&K or Hq, lines can also diagnose inhomo- 
geneities in the stellar chromosphere and photosphere, respec- 
tively ( [Bonfils et al.|2007| ). Finally, we obtained photometric ob- 
servations of a few stars to check whether plages or spots could 
produce the observed Doppler changes. 

5.1. Planetary systems 

Among the stars with clear periodic/Keplerian signals, we of 
course recover several planets that were previously known. In 
total, 14 planets are known to orbit 8 M dwarfs of our sample. 
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Nine were found by this program (Gl 176b, G1433b, G1581 b, 
c, d &e, G1667Cb and G1674b), one more was found in 1998 
by our former program on ELODIE (Gl 876 b), 2 were detected 
by concurrent programs and already confirmed by this program 
(Gl 876 c & d) and 2 were detected by concurrent programs and 
confirmed in this paper (Gl 832 b and Gl 849 b). 

• Gl 176 : From HEX radial-velocity data, it was proposed 
that it hosts a msini = 24 M® planet in 10.2-d orbit (Endl 
2QQ8| ). Howeve r, we found our HAR PS data incompat- 



et al. 



ible with that planet ( [Forveille et al.|[2Q09|) . Instead, we have 
shown evidence for a lower-mass planet with a shorter period 
(m sin / = 8 Me; P = 8 d). Much like Gl 674 (see below), Gl 176 
is a moderately active M dwarf. We also observe a second pe- 
riodic signal (P ~ 40 d), that has marginally higher power than 
the 8-d signal in our periodogram. Thanks to photometric ob- 
servations and to spectroscopic indices measured on the same 
spectra, we ident ified the 4Q-d signal as due to a spot rather than 
a second planet ( Forveille et al.|20Q9 see). We note that our sys- 



tematic Keplerian search for planets converges to solutions with 
diff'erent periods and very high eccentricities. This is mostly be- 
cause the signal associated with the 40-d period might be poorly 
described by a Keplerian function with a large eccentricity. The 
method converges to the same periods as the periodogram anal- 
ysis when we restrict the range of eccentricities to < 0.6. 



• Gl 433 : This nearby M2V dwarf is rather massive (M^ = 
0.48 Mq) for our sample. The periodogram of our HARPS RVs 
shows a clear power excess at a 7.2-d period. We failed to find a 
counter part to that signal in our activity indicators. Also, based 
on the intensity rather than the variability of Ha and Can lines, 
the star seems to have a weak magnetic activity, and most prob- 
ably a low rotational velocity. It is therefore likely that a planet 
revolves around G1433 every 7.2 days. A minimization of 
Harps RVs lead to a minimum mass of m sin / = 6 M® for that 
planet. We note that |Zechmeister et al. ( 2009| ) have reported 
UVES radial velocities for G1433 and found no significant pe- 
riodicity. The semi-amplitude of our solution 3.5+0.4 m/s trans- 
lates to an r.m.s. of 5.0+0.6 m/s that is nonetheless compatible 
with the 4.4 m/s r.m.s. reported in Z09. Also, we found Z09's 
RVs compatible with our data, provided that we use a model 
composed of a Keplerian plus a low order polynomial to fit the 
merged data sets. We refer the reader to Delfosse et al. (in prep.) 
for a detailed description. 



• Gl 581 : That system coun ts at least four p l anets, which we 
reported in three papers : ,Bonfils et al. ( 2005| ); Udry & Santos 
( |2007| ); [Mayor et al.|d2009|). We also pe rformed a s tability anal- 
ysis fo r the system in fBeust et al. ( 2008| ), updated in' Mayor et al. 
(|2009|). Composed of one Neptune-mass planet and three super- 
Earths, the system is remarkable because it includes both the first 
possibly habitable ex oplanets (G1581c &d - Selsis et aL]|2007[ 
von Bloh et al.|2007] ) and the lowest-mass planet known to date 
(G1581e - msini =1.9 M® ). Besides, we found that stability 
constrains its configurations. In the coplanar cases, inclinations 
lower than ~ 40° induce too strong interactions between "b" 
and "e", and "e" is ejected in a few hundred thousands years. 
A lower bound on the inclinations translates to upper bounds 
for planetary masses. G1581e, for instance, would not be more 
massive than ~ 3 M® if the system is coplanar. In 2010, |Vogt 



between G1581 c and d. However, we most recently demonstrate 
those planets do not exist (Forveille et al., submitted). 

• Gl 667C : This is a M2V dwarf we have intensively ob- 
served. We find several coherent signals in RV data and pos- 
sibly identify the rotation period on FWHM measurements of 
the cross-correlation function. Fitting a 1 -planet model plus a 
~ 1 .8 m s"^ yr"^ linear drift to account for the A+B stellar-binary 
companion to Gl 667C, converges toward a minimum mass of a 
super-Earth (msin/ = 5.9M®) on a short-period orbit (7.2 day). 
Adding one more planet to the model makes the fit converge to- 
ward a ~ 180 day period and a very eccentric solution, while 
the power excess identified in the periodogram of the first model 
residuals is located around 90 day. Finer analysis actually inter- 
prets that second signal as a possible harmonic of a (half-)yearly 
systematic aff'ecting few data points (Delfosse et al., in prep.). 
Filtering the signals of 2 planets plus a linear drift reveals an- 
other candidate planet (P3 = 28 d; msin / = 3.4 M®). This can- 
didate receives about 90% the amount of light received by Earth 
in our Solar System and we speculate the planet is a habitable 
canidate (see Delfosse et al., in prep., for a detailed description). 



• Gl 674 : Only 4.5 pc away from our Sun, this M2.5 dwarfs 
hosts at least one low-mass planet (m sin / = 1 1 M®; P = 4.7 d- 
Bonfils et al.|2007| ). Although a second periodic signal exists for 
Gl 674 (P2 ~ 35 d), analysis of spectroscopic indices and photo- 
metric observations shows that this additional signal originated 
from stellar surface inhomogeneities. Today and with additional 
measurements, subtracting the 4.7-d periodic signal and com- 
puting a periodogram of the residuals shows power excess at a 
period of ~ 25 d instead of 35. If due to a planet, it would corre- 
spond to a super-Earth in Gl 674 's habitable zone. However, the 
semi-amplitude K of the Keplerian orbit of that second planet 
is ~ 3.8 m/s, significantly above the residuals around the 2007 
combined fit (r.m.s.~80 cm/s). There are two diff'erent interpre- 
tations for this apparent inconsistency : either the 2007 's solu- 
tion excludes the present solution and today's 25 -d periodicity 
is spurious or, the 2007 's fit absorbed both the 35- and the 25-d 
signals, simultaneously. 



et al.| have proposed that 2 more planets orbit Gl 581, one having 
m sin / = 3.1 M® and being in the middle of the habitable zone. 



The Keplerian analysis presented in Sect. |4.4| find a lower 
significance of only ~ 94.6% for the second signal. We nev- 
ertheless apply further diagnostics. Restricting the data set to 
the 2007' RVs, we try a 1 planets sine wave model instead of 2 
Keplerians and found almost equally low residuals (r.m.s -1.1 
m/s). That strongly opposes the presence of an additional planet 
with a 3.8 m/s semi- amplitude with short or moderate orbital pe- 
riod. On the other hand, periodograms of the Yia and Can H+K 
indexes continue to peak at a ~ 35-d period (Fig. |4]), indicating 
that the signal remains coherent for these indicators. A decorre- 
lation between spectral indexes and RVs is then hard to explain. 
We conclude that the case for an additional planet is not strong 
enough with the present data set, and that data gathered after 
Apr. 1st 2009, will be necessary to conclude. 



• GI832 : A decade long RV campaign with the Anglo 
Australian Telescope (AAT) has revealed that G1832 hosts a 
long-period companion with a period almost as long as that sur- 
vey (~ 9.5 yr - [Bailey et al.|2009| ). The best fit to AAT data lead 
to a minimum mass m sin / = 0.64 Mjup. Our HARPS data do not 
span as much time and, while they do confirm with a high con- 
fidence level the long-period RV variation, they can not confirm 
the planetary nature of Gl 832b by themselves. Together with 
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Fig. 4. Periodograms of Ca II H+K (Top) and Ha (Bottom) in- 
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Fig. 5. Best solution for the model 1 planet, with A AT (blue) 
and HARPS (red) RV for Gl 832. 



the AAT dat43 our HARPS RVs refine the orbit of Gl 832b. A 
Keplerian fit using Y orbit converges to m sin / = 0.62+0.05 Mj^p 
and P = 3507 ± 181 d (Fig.jSJ. Our Keplerian analysis (§|4]fj) 
finds a possible second signal with a 35-d period. That second 
signal however, does not reach a 99% significance level in either 
our periodicity nor our Keplerian analysis. Some power excess is 
seen around ~40 d in the BIS periodogram, uncomfortably close 
to the possible 35-d periodicity. We keep measuring Gl 832 to 
clarify whether a second periodic signal is present and assess its 
true nature. 



• Gl 849 : This M3V dwarf is known since 2006 to host a 
Jupiter-mass companion ( [Butler et al.||2006[ ). The RV variation 
is clearly seen in our HARPS observations, and has no counter 



2 One RV point (with Julia n Date = 2453 243 0503) has different 
values in Table 1. and Fig. 2 of [Bailey et al.| ( [2009] ). Bailey and collab- 
orators kindly informed us of the correct value (-2.1 + 2.5 m/s) 
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Fig. 6. Best solution for the model 1 planet+drift, with Keck 
(blue) and HARPS (red) RV for Gl 849. 



part in our activity indicators (based on the shape of the cross- 
correlation function or Yia and Can spectral indices). Our obser- 
vations confirm that Gl 849 hosts a Jupiter-mass companion. 

Fitting a Keplerian orbit to HARPS observations alone con- 
verges toward a minimum mass of m sin/ = 1.17 ± 0.06 Mjup 
and a period P = 2165 ± 132 d. Together with the Keck RVs 
however, one planet is not suflftcient to explain all the RV mo- 
tion. As already suspected from Keck data, a long-term change 
is superimposed on the first periodic signal. We therefore fit 
the merged data set with a 1 planet+drift, a 2 planets and a 2 
planets+ drift model and calculate their respective FAPs. We find 
that a model more complex than 1 planet+a drift is not justified. 
For that model, our best solution ( ^|xl = 1.96) corresponds to a 
Jupiter-mass planet (m sin / = 0.99 + 0.02 Mj^p; P = 1852+19 d; 
e = 0.04 + 0.02) plus, a RV drift with a slope of -3.84 m/s/yr 
(Fig.[6| Tab.[TO]). Because Gl 849 is a nearby star (d= 8.77 + 0.16 
pc), the long-period massive companion makes it an excellent 
target for astrometric observations and direct imaging with high 
angular resolution. 



• Gl 876 : That system was known to harbor planets before our 
observations started and at that time, was even the only plan- 
etary system centered on an M dwarf. The first giant planet 
found to orbit GJ 876 was detected simultaneously by members 
of our team using the ELODIE and CORALIE spectrographs 
(f Pelfosse et ar[[T998| and by the HIRES/Keck search for exo- 
planets ( [Marcy et al.|1998| . The system was later found to host 
a second gi ant planet in a 2:1 resonance with GJ876b ( jMarcy 



et al.[[200T| ). The third planet detected around GJ876 was the 
first known super-Earth, Gl 87 6d ( (Rivera et al.[[2005| ). Because 
the 2:1 configuration of the two giant planets leads to strong 
interactions, the orbits diff'er significantly from Keplerian mo- 
tions. To model the radial velocities, one has to integrate the 
planet movement with a N-body code. This lifts the sin / de- 
generacy and measures the masses of the giant planets. A full 
N-body analysis was first performed for GJ 876 by jLaughlin &] 
[Chambers[ ( [200 1[ ) and Rivera & Lis sauer ( 2001]). From an up- 
dated set of Keck RVs, Rivera et al. ( 2005j found the third planet 
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Table 10. Fitted orbital solutions for the Gl 849 (Iplanet + a linear drift) and Gl 832 (1 planet). 



Parameter 




G J 849 b 


G J 832 b 


P 


[days] 


1845 + 15 


3507 + 181 


T 


[JD-2400000] 


54000 + 150 


54224 + 302 


e 




0.05 + 0.03 


0.08 + 0.05 


to 


[deg] 


298 + 29 


254 + 35 


K 


[m s-i] 


24.4 + 0.7 


14.1 + 1.1 


V(Keck/AAT) 


[km s"^] 


-0.0192 + 0.0012 


0.0033 + 0.0013 


ViHarps) 


[km s-i] 


-15.0896 + 0.0005 


13.3471 + 0.0013 


Slope 


[m s~^ yr~^] 


-4.76+0.33 




/(m) 


[lO-^Mo] 


2.77 


1.00 


m2 sin / 


[Me] 


0.91 


0.62 


a 


[AU] 


2.32 


3.46 


o- (O-C) (AAT/Keck) 


[ms-^] 


3.64 


4.70 


a (O-C) (Harps) 


[ms-i] 


2.08 


1.77 






1.83 


2.36 



only because planet-planet interactions were properly accounted 
in the fitting procedure. Those authors sti ll had to assume c opla- 
nar orbits to assign a mass to each planet. Bean et al. (2009) then 
combined the Keck RVs with HST astrometry to both measure 
the masses in the coplanar case and to measure the relative incli- 
nation between planets "b" and "c". IMost recently, we used our 
HARPS data and the published Keck measurements to model the 
system and measure the relative inclina tion of both giant pl anets 
(< 1°), relying on RVs only. The paper (' Correia et al.|201 ) also 
analyzed the dynamical stability and show that the librations am- 
plitude are smaller than 2° thanks to a damping process acting 
during the planet formation. 

5.2. Activity dominated variations 

We group in this section the "active" cases. They are stars tested 
positively for periodicity and/or Keplerian signal, and their mea- 
surement variability correlates with an activity indicator. We do 
not show all diagnostics for each star, but rather pick the most 
illustrative. A statistical discussion of all activity indicators will 
be presented in a separate paper (Bonfils et al. 2010, in prep.). 

• GI205 : A periodogram of the velocities identifies excess 
power around 32.8 d. Our Keplerian search with Yorbit, for a 
single planet or for a 1 planet+drift model, converges toward 
either a similar period or 0.970 d, an alias with the 1-d sam- 
pling. We find indications that the variation is intrinsic to the 
star from both spectral indices and photometric observations. 
Considering the whole dataset, we identify excess power around 
33 d in periodograms of Ha and Can H-fK indices, though, 
those peaks are not the highest. Restricting instead the dataset 
to one observational season, a strong power excess around the 
period 33 d dominates the periodogram (Fig.jTj). With addition- 
ally note a high correlation between Ha and Can H-hK indices 
(their Pearson correlation coefficient is 0.97), suggesting both 
variations originate from the same surface homogeneity. Also, 
photometric monitoring of Gl 205 reports a similar period for the 
stellar rotation (33.61 d - |Kiraga & Stepien|20Q7| . The observed 
RV modulation is most probably due to surface inhomogeneities, 
which remains stable over one season but not over several. 



• Gl 358 : We have gathered 28 measurements for Gl 358. They 
show significant variability with a periodicity of ~ 26 d. In the 
RV time series periodogram, we also identify power excess at 
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Fig. 7. For G1205, the top 2 panels show the periodograms for 
Ca II H-hK index and Ha index, including all data. The bot- 
tom 2 panels show the periodograms for Ca II H-fK index and 
Ha index, restricting the dataset to one observational season 
(2, 453, 600 < BJD < 2, 543, 900). 



the first harmonic of that period (~ 13d). The RV modulation 
is well described by a Keplerian orbit. However, we also ob- 
serve similar variability in the FWHM of the CCF as well as a 
possible anti-correlation between RV and spectral lines asymme- 
try (see Fig [8]), as measured by the CCF bisector span ( Qu eloz] 
let al.|2001| ). The Pearson's correlation for BIS and RV is -0.40, 
and rises to -0.67 for the 2007 measurement subset. Photometric 
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Fig. 8. Top : Periodogram for the FWHM of CCFs for Gl 358. 
Bottom : Possible correlation between bisector spans and RV 
data for G1358. 
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Fig. 9. Top : Periodogram for the CCFs bisector span for Gl 388. 
Bottom : Strong correlation between bisector spans and RV data 
for G1388. 



monitoring by|Kiraga & Stepien ( |2QQ7 ) found a rotational period 
of 25.26 d for Gl 358. All this favors a scenario where a stellar 
surface inhomogeneity such as a spot or a plage produces the RV 
change, rather than a planet. 

• GI388 (AD Leo) : The periodogram of its RV time series 
shows an important power excess at short-period, with two 
prominent peaks at 1.8 and 2.22 d, consistent with 2.24-d ro- 
tational period reported by |Morin et al.| ( |2008| ). These are 1- 
d aliases of each other, and the later is slightly stronger. Here, 
the bisector span demonstrates that stellar activity is responsible 
for the variation. Its periodogram shows a broad power excess 
at short period, and it is strongly anti-correlated to RV (with a 
Person's correlation coefficient of -0.81 - see Fig. [9]). Correcting 
from the BIS-RV correlation by subtracting a linear fit does de- 
crease the r.m.s from 24 to 14 m/s, but leaves some power excess 
around ~ 2 d. 



• Gl 479 We observe significant power excesses in the RV time 
series at two periods, -11.3 and 23 - 24 d, with the shorter pe- 
riod roughly half the longer one. The RVs vary with an ampli- 
tude of ~27 m/s and an r.m.s. of 4.13 m/s. Modeling that RV 
variability with Keplerians converges toward 2 planets with very 
similar periods (23.23 and 23.40 d), which would clearly be an 
unstable system. Gl 479 shares its M3 spectral type with Gl 674 
and G1581, which we use as benchmarks for their moderate and 
weak magnetic activity, respectively. From a spectral index built 
on the Can H&K lines, we find that G1479 has a magnetic ac- 
tivity intermediate to Gl 674 and Gl 581 . Neither the bisector nor 
the spectral indices show any significant periodicity or correla- 
tion with RVs. However, we have complemented our diagnos- 
tic for stellar activity with a photometric campaign with Euler 
Telescope (La Silla). The periodogram of the photometry shows 
a maximum power excess for the period 23.75 d, similar to the 
RV periodicity. The photometry phase folded to the 23.75-d pe- 
riod varies with a peak-to-peak amplitude of 5% and complex 
patterns. We cannot ascertain the rotational period of Gl 479, but 
a 5% variability can explain the observed RV variations down to 
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Fig. 10. Top : Periodogram of G1479 photometry. Bottom 
Phase-folded to the 23.75-d period. 



very slow rotation. We cannot correlate the photometry with RV 
variation because they were not taken at the same time but, phas- 
ing both to a 23 d period, we found a 0.24-phase shift consistent 
with a spot. The observed RV variability is therefore probably 
due to Gl 479 magnetic activity rather than to planets. 



• Gl 526 We observe RV periodicity with power excess close 
to 50 d and a FAP approaching 1%. As expected for a moderate 
or long period, we do find corresponding changes neither in BIS 
nor FWHM. As for Gl 674, spectral indices or photometry are 
then more informative. For Gl 526, we do find that RV is weakly 
correlated to Ha and that Ca II H-hK index varies with a clear 



50-d period (Fig. 11). Because the observed period is similar for 
the calcium index and the RV shift, we interpret that RV changes 
as due to magnetic activity rather than a planet. 
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Fig. 11. Top : Periodogram for the Ca II H+K index for Gl 526. 
Bottom : Correlation between Ha index and RV data for Gl 526. 
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Fig. 12. Correlation between bisector inverse slope and RV data 
for G1846. 



• Gl 846 : We observe RV variability with significant power ex- 
cesses in the periodogram at several periods (7.4, 7.9 and 10.6 
d), plus their aliases with the 1-d sampling, near 1-d. The bisec- 
tor inverse slope is well correlated with the RV (Fig. [12]). Like 
Gl 358 and Gl 388, Gl 846 is also a clear case of stellar intrinsic 
variability rather than planetary-companion Doppler shift. 



5.3. Unclear cases 

• Gl 273 This M3.5 dwarf shows RV variability and tests pos- 
itively for a RV drift according to the -probability test (and 
approaches FAP=1% for the permutation test). The RV pe- 
riodogram indicates significant power excess at a period of 
~ 434 d, and Yorbit find a good solution for a planet with 
P ~ 440 d, with or without a supplementary drift. However, 
good solutions are found with uncomfortably high eccentricity 



and, most importantly, a poor phase coverage. Among activity 
indicators, only the Ha index has power excess for long peri- 
ods but with a diff'erent period (~ 500 - 600 d). Once a linear 
drift is subtracted, 2 RV points stand out (with BJD=2,454,775 
and BJD=2,454,779). They have a value 8-10 m/s lower than 
the residual mean. If we fit a drift again on the original data, 
considering all but those 2 points, the 440-d power excess dis- 
appear from the periodogram of the residuals. This suggest that, 
if these 2 particular points were outliers, the 440-d period could 
be an alias between a long-term RV drift and the 1-yr sampling. 
However, a direct inspection of Gl 273 spectra, cross-correlation 
functions, and spectral indexes give no reason to exclude those 
values. We conclude that a firm conclusion on G1273 would be 
premature and will obtain more data. 

• Gl 887 Formally significant models are found with 1 and 3 
planets but all have converged to solutions with very high or 
unrealistic eccentricities. Most probably, the RV variability of 
Gl 887 does not match a Keplerian motion, and our automatic 
search got confused. Also, we do not find significant periodicity 
or correlation with RVs in our diagnostics. 



6. Detection limits 

While the previous sections focuses on signal detection, the 
present one aims at giving upper limits on the signals we were 
not able to detect. For individual stars, the upper limit translates 
into which planet, as a function of its mass and period, can be 
ruled out given our observations. For the sample, all upper limits 
taken together convert into a survey eflftciency and can be used 
to measure statistical properties (Sect. [7]). 

To derive a period-mass limit above which we can rule out 
the presence of a planet, given our observations, we start with the 
periodogram analysis presented in |4.3| We make the hypothesis 
the time series consists of noise only. Like for the FAP calcula- 
tions, we evaluate the noise in the periodogram by generating 
virtual data sets. The virtual data are created by shuffling the 
time series while retaining the observing dates (i.e. by bootstrap 
randomization). And for each trial we compute a periodogram 
again. This time however, we do not look for the most powerful 
peak. We rather keep all periodograms and build a distribution 
of powers, at each period. For a given period then, the power 
distribution tells us the range of power values compatible with 
no planet, i.e. compatible with our hypothesis that the time se- 
ries consists of noise only. In the same manner as for the FAP 
computation, we can also evaluate the probability that a given 
power value occurred by chance just by counting the fraction of 
the power distribution with lower values. 

Once we know the power distribution compatible with no 
planet, we inject trial planetary orbits into the data. In this pa- 
per, we restrict our analysis to circular orbits and therefore fix 
the eccentricity and argument of periastron to zero. We never- 
theless note that eccentricities as hi gh as 0.5 do not aflfect much 
the upper limit on planet detection ( Endl et al.|2002[ Cunmiing 
[& Dragomir 2009). We thus add sine waves, choosing a period 

a semi-amplitude K and a phase T. We explore all periods 
computed in the periodograms, from 1.5 to 10,000 day, with a 
linear sampling in frequency of step 1/20,000 day"\ and for 12 
equi- spaced phases. At each trial period, we start with the semi- 
amplitude of the best sine fit to the original time series Kobs^ and 
compute the new periodogram power psim for that period. We 
next increase the semi-amplitude K until Psim reaches a value 
with a probability as low as or lower than 1%, if the data were 
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noise only. On one hand, we impose that power threshold for 
all our 12 trial phases and obtain a conservative detection limit. 
On the other hand, we averaged that power threshold over our 
12 trial phases and obtain a statistical (or phase -averaged) de- 
tection limit. Eventually, for both detection limits, we convert 
the K semi-amplitude to planetary mas^ and orbital period to 
orbital semi-major axij^ using the stellar mass listed in Tablejs] 

The method d escribed above ha s been applied before us by 
Gumming et~aLl ([1999, 2008) and [Zechmeister et al.| ( [20091 ). 
We note a small difference between Gumming et al's and 
Zechmeister et al's approaches. On one hand. Gumming and col- 
laborators add the trial sine wave to normally distributed noise 
and choose as variance the r.m.s around the best sine wave to 
the observed data. On the other hand, Zechmeister and collabo- 
rators choose not to make trial versions for the noise and consid- 
ered the observed data as the noise itself (to which they add the 
trial orbit). We choose Zechmeister et al.'s approach because we 
believe it is more conservative when, in some cases, the noise 
departs from a normal distribution. 

We report both conservative and phase-averaged upper lim- 
its for the 98 time series with more than 4 measurements (an 



example is shown in Fig . [13] an d [T4| for Gl 581 and we group the 
figures of all stars in Fig. |18| and [19] which is only available 
online). When a periodic variation has been attributed to stel- 
lar magnetic activity, we know the main variability is not due to 
a planet. We therefore apply a first order correction to the RV 
data by subtracting the best sine fit. We choose a simple sine 
wave rather than a more complex function (like a Keplerian) 
because the fundamental Fourier term is the least informative 
choice, and therefore the most conservative. When instead, we 
have identified the RV variability is due to one or more plan- 
ets, we are interested in the upper-limit imposed by the residuals 
around the solution. We therefore subtract the best Keplerian so- 
lution to the time series before computing the upper limit. For the 
multi-planetary systems G1581, G1667G and G1876, we com- 
pute the periodogram once with raw time series and once with 
the residuals around the full orbital solution (with all detected 
planets). Because the giant planets in Gl 876 system undergo 
strong mutual interactions, we use a N-body integration to com- 
pute the residuals ( [Gorreia et al.||2010p . For G1674 and Gl 176 
(that show both planet- and activity-induced variation), we use 
a 1 Keplerian sine model to fit the RVs. At last, when we ob- 
served no variability, variability without periodicity, or periodic 
variability without a well identified cause (planet v^. magnetic 
activity), we use the raw time series to compute periodograms 
and upper limits. 

7, Planet occurrence 

All together, the phase- averaged detection limits calculated for 
individual stars give the survey efficiency, which is eventually 
used to correct for the detection incompleteness and derive pre- 
cise occurrence of planets. Although the statistical analysis of 
our survey is the purpose of a companion paper (Bonfils et al., 
in prep.) we give here a first aper^u. 

In an m sin /-period diagram, we pool together the phase- 
averaged detection limits computed in Sect. [6] and, for each pe- 
riod, count the number of ruled out planets more massive than 
a given mass (with our 99% criterion). This synthesis is shown 
in Fig. [T5] together with iso-contours for 1, 10, 20, 30, 40, 50, 
60, 70, 80, and 90 stars. We also overlay the planet detections 



including all planets of multi-planet systems. Such a diagram 
is especially useful to compare the planet occurrence for dif- 
ferent domains of masses and periods. For instance, for peri- 
ods P < 100 day, our sample counts 1 host star (G1876) with 
two giant planets (msin/ = 0.5-10 Mjup) but 7 super-Earths 
(msin/ = 1-10 M®). Whereas, our survey is sensitive to 
Gl 876b-like planets for 92 stars, short-period super Earths could 
be detected for only -5-20 stars. It is therefore obvious that 
super-Earths are much more frequent than Gl 876b-like giants. 

For more precise estimates we delineate regions in the mass- 
period diagram and approximate the frequency of planets by the 
ratio / = NdlN^^eff^ where Nd is the number of planets detected 
in that region and where N^^^eff is the number of stars whose 
detection limits confidently exclude planets with similar mass 
and period. We evaluate A^*,^// with Monte-Garlo sampling : 
we draw random mass and period within the region delineated 
(assuming a log-uniform probability for both quantities), use 
the m sin /-period diagram of Fig. [15] to give a local estimate of 
Ni^^eff and, with many trials, compute an averaged N^^^eff value. 
Table. 1 1 reports the number of detections, averaged values for 
N^^eff and corresponding occurrence of planets for different re- 
gions chosen in Fig. [15 



Numbers of Table. 



1 1 do confirm that planets are increas- 



a = {PI27ifl\GM^yi^ 



ingly abundant toward lower-mass and longer-period planets. 

At last, we estimate 77®, the frequency of habitable planets 
orbiting M dwarfs. For each star, we use the locations for both 
the inner (aHzjn) and outer (aHz,out) edges of the habitable zone 
computed in Sect. [2] and we consider that habitable planets have 
msin / between 1 and 10 M®. To evaluate the sensitivity of our 
survey to habitable planets, we compute a new A^*,e// for the 
habitable zone with a Monte-Garlo approach again. We draw 
random masses between 1-10 M® and random semi-major axis 
between anzjn and aHz,out^ choosing a log-uniform probability 
for both mass and semi-major axis. We screened the detection 
limits computed previously and increment A^*,^// when the trial 
fall above that threshold. After normalizing A^*,e// by the num- 
ber of trials we found A^*,e// = 4.84. As among our detections 
two planets (G1581d and G1667Gc) falls in the habitable zone, 
we have Nd = 2 and therefore 77® = 0.41^q^3. 

Alternatively, we measure that 1 1 (resp. 3) stars of our sam- 
ple have time-series precise enough to detect planet with same 
mass and period as G1581d (resp. G1667Gc), which lead to a 
very similar estimate of 77® (~ 42%). 

8. Conclusion 

We have reported on HARPS guaranteed-time observations for a 
volume-limited sample of nearby M dwarfs. The paper develops 
on a systematic analysis of the time series for 102 M dwarfs. 
It analyzes their variability, look for possible trends, and search 
them for periodicities and Keplerian signals. 

We find significant periodic signals for 14 stars and linear 
trends for 15. We recover the signal for 14 known planets in 8 
systems. In particular, we confirm the detection of 2 giant plan- 
ets, G1849b and G1832b, and we conffim that an additional 
long-period companion is probably orbiting Gl 849. We analyze 
the RV periodicity against stellar diagnostics for 8 other stars, 
and find evidences the observed RV variation originate from stel- 
lar surface inhomogeneities for all but one (Gl 273). We find pe- 
riodic RV variation in Gl 273 time series without counter part in 
activity indicators, though the phase coverage is too poor for a 
robust detection. 

Our search for planets with HARPS has detected 9 planets 
in that sample alone, and a total of 1 1 planets when counting 2 



14 X. Bonfils et al.: The HARPS M-dwarf sample 

Table 11. Occurrence of planets around M dwarfs for various regions of the m sin /-period diagram. 
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[day] 
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Neff = 13.97 


Neff = 5.79 


Neff = 1.53 


Neff = 0.03 






/ = 0.52^^0 



M dwarfs from another complementary sample. Our detections 
includes the lowest-mass planet known so far and the first proto- 
type of habitable super-Earths. They are the fruit of slightly less 
than 500 h of observing time on a 3.6-m telescope, now-day s 
considered as a modest size telescope. 

Beyond individual detections, we have also reported a first 
statistical analysis of our survey. We derived the occurrence of 
M-dwarf planets for diff'erent regions of the m sin /-period dia- 
gram. In particular, we find that giant planets (msin/ = 100 - 
1, 000 Me) have a low frequency (e.g. / < 1% for P = 1 - 10 d 
and / = 0.02+^;^^ for P = 10 - 100 d), whereas super-Earths 
(m sin / = 1-10 M®) are likely very abundant ( / = 0.36^^^^ for 
P = 1 - 10 d and / = 0.35+^ for P = 10 - 100 d). We also ob- 
tained 7/® = 0.41^Q^3, the frequency of habitable planets orbiting 
M dwarfs. Considering M dwarfs dominate the stellar count in 
the Galaxy, that estimate weighs considerably in the measure of 
the frequency of habitable planets in our Galaxy. Also, for the 
first time, 7/® is a direct measure and not a number extrapolated 
from the statistic of more massive planets. 

Of course, much refinements are possible. For instance, back 
to Fig. [T] we indicate with vertical ticks (above the histograms) 
the V and mass values for the known planet-host stars included 
in our sample. It is striking that all planet-host stars are found in 
the brightest and more massive halves of the two distributions. 
This is reminiscent of what we observe between solar type stars 
and M dwarfs. Solar-like stars have many more detected plan- 
ets, but observational advantages and disadvantages compared 
to M dwarfs are difficult to weight. On the one hand solar-type 
stars a re brighter and have a lower jitter level (e.g. Hartman et al. 
201 1| ), but on the other hand they are more massive and the reflex 
motion induced by a given planet is weaker. To know whether we 
face an observational bias or a true stellar mass dependance in 
the formation of planets, we need to evaluate the detection effi- 
ciency for all mass ranges, which is the purpose of a forthcoming 
companion paper (Bonfils et al., in prep.). 

Acknowledgements. We express our gratitude to Martin Kuerster who serve as 
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the manuscript. 



References 

Apps, K., Clubb, K. I., Fischer, D. A., et al. 2010, PASP, 122, 156 
Bailey, J., Butler, R. P., Tinney, C. G., et al. 2009, ApJ, 690, 743 
Baranne, A., Queloz, D., Mayor, M., et al. 1996, A&AS, 119, 373 
Batista, V., Gould, A., Dieters, S., et al. 2011, Astronomy & Astrophysics, 529, 
102 



Bean, J. L., Seifahrt, A., Hartman, H., et al. 2009, eprint arXiv, 0912, 3, submit- 
ted to ApJ Letters 
Beaulieu, J.-R, Bennett, D. R, Fouque, R, et al. 2006, Nature, 439, 437 
Benedict, G. R, McArthur, B. E., Forveille, T., et al. 2002, ApJ, 581, LI 15 
Bennett, D. R 2009, arXiv, astro-ph.EP 

Bennett, D. R, Bond, I. A., Udalski, A., et al. 2008, ApJ, 684, 663 
Beust, H., Bonfils, X., Delfosse, X., & Udry, S. 2008, A&A, 479, 277 
Boisse, I., Bouchy, R, Chazelas, B., et al. 2010, eprint arXiv, 1001, 794 
Bond, I. A., Udalski, A., Jaroszyhski, M., et al. 2004, ApJ, 606, L155 
Bonfils, X., Delfosse, X., Forveille, T., Mayor, M., & Udry, S. 2007, Proceedings 

of the conference In the Spirit of Bernard Lyot: The Direct Detection of 

Planets and Circumstellar Disks in the 21st Century. June 04 - 08, 21 
Bonfils, X., Delfosse, X., Udry, S., Forveille, T., & Naef, D. 2006, Tenth 

Anniversary of 51 Peg-b: Status of and prospects for hot Jupiter studies. 

Colloquium held at Observatoire de Haute Provence, 111, iSBN: 2-914601- 

174 

Bonfils, X., Delfosse, X., Udry, S., et al. 2005, A&A, 442, 635 
Bonfils, X., Gillon, M., Forveille, T, et al. 2011, Astronomy & Astrophysics, 
528, 111 

Bouchy, R, Pepe, R, & Queloz, D. 2001, A&A, 374, 733 
Butler, R. R, Johnson, J. A., Marcy, G. W., et al. 2006, PASP, 118, 1685 
Butler, R. R, Vogt, S. S., Marcy, G. W., et al. 2004, ApJ, 617, 580 
Chabrier, G. & Baraffe, I. 2000, ARA&A, 38, 337 
Charbonneau, D., Berta, Z. K., Irwin, J., et al. 2009, Nature, 462, 891 
Charbonneau, D., Brown, T. M., Latham, D. W., & Mayor, M. 2000, ApJ, 529, 
L45 

Chauvin, G., Lagrange, A.-M., Dumas, C, et al. 2004, A&A, 425, L29 
Correia, A. C. M., Couetdic, J., Laskar, J., et al. 2010, A&A, 511, 21 
Gumming, A., Butler, R. R, Marcy, G. W., et al. 2008, PASP, 120, 531 
Gumming, A. & Dragomir, D. 2009, eprint arXiv, 0909, 4758, 16 pages, 9 fig- 
ures, to appear in MNRAS 
Gumming, A., Marcy, G. W., & Butler, R. R 1999, ApJ, 526, 890 
Cutri, R. M., Skrutskie, M. R, van Dyk, S., et al. 2003, The IRSA 2MASS All- 
Sky Point Source Catalog 
Delfosse, X., Forveille, T, Beuzit, J.-L., et al. 1999, A&A, 344, 897 
Delfosse, X., Forveille, T, Perrier, C, & Mayor, M. 1998, A&A, 331, 581 
Delfosse, X., Forveille, T, Segransan, D., et al. 2000, A&A, 364, 217 
Deming, D., Harrington, J., Laughlin, G., et al. 2007, ApJ, 667, L199 
Demory, B.-O., Gillon, M., Barman, T, et al. 2007, A&A, 475, 1125 
Dong, S., Bond, I. A., Gould, A., et al. 2009a, ApJ, 698, 1826 
Dong, S., Gould, A., Udalski, A., et al. 2009b, ApJ, 695, 970 
Endl, M., Cochran, W. D., Kurster, M., et al. 2006, ApJ, 649, 436 
Endl, M., Cochran, W. D., Wittenmyer, R. A., & Boss, A. R 2008, ApJ, 673, 
1165 

Endl, M. & Kurster, M. 2008, A&A, 488, 1 149 
Endl, M., Kurster, M., Els, S., et al. 2002, A&A, 392, 671 
Forveille, T, Bonfils, X., Curto, G. L., et al. 2011, A&A, 526, 141 
Forveille, T, Bonfils, X., Delfosse, X., et al. 2009, A&A, 493, 645 
Gatewood, G. & Eichhorn, H. 1973, Astron. J., 78, 769 
Gatewood, G. D. 1995, Ap&SS, 223, 91 

Gaudi, B. S., Bennett, D. R, Udalski, A., et al. 2008, Science, 319, 927 
Gillon, M., Demory, B.-O., Barman, T, et al. 2007a, A&A, 471, L51 
Gillon, M., Pont, R, Demory, B.-O., et al. 2007b, A&A, 472, L13 
Gliese, W. & JahreiB, H. 1991, On: The Astronomical Data Center CD-ROM: 

Selected Astronomical Catalogs, aRIPRINT Id. #1991016 
Gould, A., Udalski, A., An, D., et al. 2006, ApJ, 644, L37 



X. Bonfils et al.: The HARPS M-dwarf sample 



15 



Haghighipour, N., Vogt, S. S., Butler, R. R, et al. 2010, ApJ, 715, 271 
Halbwachs, J. L., Arenou, R, Mayor, M., Udry, S., & Queloz, D. 2000, A&A, 
355,581 

Hartman, J. D., Bakos, G. A., Noyes, R. W., et al. 2011, AJ, 141, 166 

Henry, G. W., Marcy, G. W., Butler, R. P., & Vogt, S. S. 2000, ApJ, 529, L41 

Henry, T. J., Jao, W.-C, Subasavage, J. P, et al. 2006, AJ, 132, 2360 

Hershey, J. L. 1973, Astron. J., 78, 421 

Ida, S. & Lin, D. N. C. 2005, ApJ, 626, 1045 

Janczak, J., Fukui, A., Dong, S., et al. 2010, ApJ, 711, 731 

Joergens, V. & MuUer, A. 2007, ApJ, 666, LI 13 

Johnson, J. A., Apps, K., Gazak, J. Z., et al. 2010, eprint arXiv, 1008, 4141 
Johnson, J. A., Butler, R. P, Marcy, G. W., et al. 2007, ApJ, 670, 833 
Kalas, P, Liu, M. C., & Matthews, B. C. 2004, Science, 303, 1990 
Kasting, J. E, Whitmire, D. P, & Reynolds, R. T. 1993, Icarus, 101, 108 
Kiraga, M. & Stepien, K. 2007, Acta Astronomica, 57, 149 
Kubas, D., Beauheu, J. P, Bennett, D. P, et al. 2010, eprint arXiv, 1009, 5665 
Kiirster, M., Endl, M., Rouesnel, R, et al. 2003, A&A, 403, 1077 
Laughlin, G. & Chambers, J. E. 2001, ApJ, 551, L109 

Leggett, S. K., AUard, E, Geballe, T. R., Hauschildt, P H., & Schweitzer, A. 

2001, ApJ, 548, 908 
Lomb, N. R. 1976, Ap&SS, 39, 447 

Maness, H. L., Marcy, G. W., Ford, E. B., et al. 2007, PASP, 119, 90 
Marcy, G. W., Butler, R. P, Fischer, D., et al. 2001, ApJ, 556, 296 
Marcy, G. W., Butler, R. P, Vogt, S. S., Fischer, D., & Lissauer, J. J. 1998, ApJ, 
505, L147 

Mayor, M., Bonfils, X., Forveille, T., et al. 2009, A&A, 507, 487 

Mayor, M., Pepe, E, Queloz, D., et al. 2003, The Messenger, 114, 20 

Mayor, M. & Queloz, D. 1995, Nature, 378, 355 

Mazeh, T, Naef, D., Torres, G., et al. 2000, ApJ, 532, L55 

Miyake, N., Sumi, T, Dong, S., et al. 2010, eprint arXiv, 1010, 1809 

Montagnier, G., Segransan, D., Beuzit, J.-L., et al. 2006, A&A, 460, L19 

Morin, J., Donati, J.-E, Petit, P, et al. 2008, MNRAS, 390, 567 

Murdoch, K. A., Hearnshaw, J. B., & Clark, M. 1993, ApJ, 413, 349 

Nutzman, P & Charbonneau, D. 2008, PASP, 120, 317 

Pepe, E, Mayor, M., Galland, E, et al. 2002, A&A, 388, 632 

Pepe, E, Mayor, M., Queloz, D., et al. 2004, A&A, 423, 385 

Queloz, D., Henry, G. W., Sivan, J. P, et al. 2001, A&A, 379, 279 

Reid, I. N., Hawley, S. L., & Gizis, J. E. 1995, AJ, 110, 1838 

Rivera, E. J. & Lissauer, J. J. 2001, ApJ, 558, 392 

Rivera, E. J., Lissauer, J. J., Butler, R. P, et al. 2005, ApJ, 634, 625 

Scargle, J. D. 1982, ApJ, 263, 835 

Schlesinger, E 1917, AJ, 30, 137 

Selsis, E, Kasting, J. E, Levrard, B., et al. 2007, A&A, 476, 1373 

Smart, R. L., loannidis, G., Jones, H. R. A., Bucciarelli, B., & Lattanzi, M. G. 

2010, A&A, 514, 84 
Sumi, T, Bennett, D. P, Bond, I. A., et al. 2010, ApJ, 710, 1641 
Udalski, A., Jaroszyhski, M., Paczyhski, B., et al. 2005, ApJ, 628, L109 
Udry, S. & Santos, N. C. 2007, Annual Review of Astronomy & Astrophysics, 

45, 397 

van Altena, W. E, Lee, J. T, & Hoffieit, E. D. 1995, New Haven 
van de Kamp, P 1963, Astron. J., 68, 515 

van de Kamp, P 1982, Vistas in Astronomy (ISSN 0083-6656), 26, 141 
van Leeuwen, E 2007, A&A, 474, 653 

Vogt, S. S., Butler, R. P, Rivera, E. J., et al. 2010, eprint arXiv, 1009, 5733 
von Bloh, W., Bounama, C, Cuntz, M., & Franck, S. 2007, A&A, 476, 1365 
Walker, G. A. H., Walker, A. R., Irwin, A. W., et al. 1995, Icarus, 116, 359 
Zechmeister, M. & Kurster, M. 2009, A&A, 496, 577 
Zechmeister, M., Kurster, M., & Endl, M. 2009, A&A, 505, 859 



16 



X. Bonfils et al.: The HARPS M-dwarf sample 



I 

S -10 







1 1 




ml 

m « 




W. 














1 1 


1 


1 1 1 

GI54.1 


; 

1...?.. 


1 s 












< 








ro 
op 

U 




1 


f . 
















H 
VO 

.. o. . . 

H 

5* 


£ I 


I 

I. 




















GI176 


if 


3— 

iff 














1 1 


1 


! 1 1 

GI203 


T 

I ^ 


J I. 












.. m . 

(N 
H 

■ CO 




1 1 


1 




1 - 




- I 












1 1 


I 1 1 

GI273 




. v# i 






X 


















vo 
vo 

.. Q . . 
(N 








1 


1 1 


1 1 


1 



1 


1 1 


1 1 


1 


(N 

o 

.. o. . . 

H 


I' 


1 


i- 










1 


1 1 


1 1 


1 


1 1 1 

L707-74 


J 




h 














- p4 








00 

"!t 

. . iH - - 
(A 
Z 

- - J - - 




[ I 


I 


.1. 

i.. 
















in 
vo 

.. Q . . 








H 

5* 


I 


1 












- -fn 

(N 

tN 

.. .r4 - - 
(/) 

z 

- - J - - 












1 








1 




1 


I 1 1 

GI205 












1 1 1 

1IP3:1292 


1 1 


1 1 


1 


I 














1 


1 1 


1 1 


1 


CO 
Oi 
. . iH - - 
(A 
Z 

- - J - - 




: i 


i 


i- 
















! 1 1 

GJ1:123 


1 


, I 

' '. 


l-l- 

[ 




1 


1 i 1 


1 



1 


1 1 


1 1 


1 


1 — I — 

GI:12 




I 


i- 










1 


1 1 


1 1 


1 


1 1 1 

Gl:87 


if 






»- - 












- -m 








H 

in 

. . iH - - 
(A 
Z 

- - J - - 


.n. 


i 


















00 

vo 

.. o. . . 

H 




I 


' 1 














- p4 
CO 
IV 

- -r4 - 

(A 

z 

- - J - - 














i ^ 


X 


! 








1 


1 1 


1 1 


1 


I 1 1 

GI213 






T 

X 
S. 










1 

- I-+ 

(N 

.. .00. . . 

o 

H 


1 1 


1 1 


1 


I ^ 






i 












1 


1 1 


1 1 


1 


I 1 1 

GI299 


u 






I 


I 




















! 1 1 

GI341 




3 

1 




E 


1 


1 1 


1 1 


1 



- «!t 
CO 
H 

. . .iH . . 
(/) 

z 

.. . J . . 












iJ. 


T 

i 




1 


1 1 


1 1 


1 


1 1 1 

GIlbSB 












i 


















! 1 1 

GI145 


I " 


1 





















1 




66C 


1 






H 




X 



















H 








- r4 






















I 1 1 

GI229 
























« . . . _ 


OQ 

o 

.. in 


f i 




CM 










i 







500 



1500 



2500 



500 



1500 



2500 



500 



1500 



2500 



500 



5 ii 



1500 



i 



2500 



BJD-2,452,000.0 [day] 



Fig. 3. Radial- velocity time series. 



X. Bonfils et al.: The HARPS M-dwarf sample 



17 



> 







1 1 


1 


in 

- £2 


'l 

i 


1 " 


z 













1 1 1 

« 


1 


CM 
00 

- rn 
13 












z 2 

i 












vo 
o 
- ^ 


t 


I 






I 

-L 





















in 
vo 
- ^ 




i - 




Ir 


















vo 

00 

- 




















1 1 


1 

*l 


H 

in 
- ift 


• t 






i--- 










1 1 


1 

I 

z 


1 1 

GI618A 


3 

Z ] 


E 




i 






*. . . 


1 1 


1 


- VO 

iD 


; 




M - - - - 


















vo 

00 

- vo 
iD 








1 


1 1 


1 1 


1 



1 


1 1 


1 1 


1 


1 1 

GI358 


i - 


1 












1 


1 1 


1 1 


1 


1 1 1 

GI393 






















! 1 1 

GI433 


















! 1 1 

: GI479 : 




Hi 
















H 

- m 

















1 


1 1 


1 1 


1 


! 1 1 

: GI555 
















1 


1 1 


1 1 


1 


1 1 

GI628 


z 

* a 


■ *l 


f 






i 


1 


1 1 


1 1 


1 


! ! 1 

GI67:8.1/ 






















! 1 1 

GI693 




■ -5 


_z_ _ _ 


1 


1 1 


1 1 


1 



1 


1 1 


1 1 


1 


1 1 

GI367 
















1 


1 1 


1 1 


1 


1 1 1 

LHS288 








( 


1 
























! 1 1 

GI438 




























! 1 1 

LHS337 


1... 






















GI526 

















1 


1 1 1 


1 


! 1 1 

GI569A 




I 












1 


1 1 1 


1 


1 1 

GI643 


I- 














1 


1 1 


1 1 


! 1 1 

GI680 










i 








! 1 1 

GI699 






1 


1 1 


1 i 











(N 
H 
" H 

5* 




i 












1 


1 1 


1 1 


1 


GI402 


-i -2- 








J. 




















! 1 1 

GI447 










Z * * 




















! 1 1 

GI480.1 


I 





















GI536 


z 


|z 









1 


1 1 


1 1 


1 


I 1 1 

GI588 


- ! 


?« - 













1 


1 1 1 


u 








.. vo. . . 




vo 













1 


1 1 


1 1 


1 


I 1 1 

GI682 




z 




S 



















! 1 1 

GI701 


"at " " : 

z 


.. * 




1 


1 1 


1 1 


1 



500 



1500 



2500 



500 



1500 



2500 



500 



1500 



2500 



500 



1500 



2500 



BJD-2,452,000.0 [day] 



Fig. 3. Radial- velocity time series (continued). 
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Fig. 3. Radial- velocity time series (continued). 
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Fig. 13. Conservative detection limit applied to G1581. Planets with minimum mass above the limit are excluded with a 99% 
confidence level for all 12 trial phases. The upper curve shows the limit before any planetary signal is remove to the RV time series. 
The sharp decrease in detection sensitivity around the period 5.3 day is caused by the RV signal of Gl 581b. The lower curve shows 
the limit after the best 4-planet Keplerian fit has been subtracted. The sharp decrease in sensitivity around the period 2 day is due to 
sampling. Venus and Mars criterion delineate the habitable zone, shown in blue. The vertical yellow dashed line marks the duration 
of the survey. 
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Fig. 14. Phase-averaged detection limit applied to G1581. Planets with minimum mass above the limit are excluded with a 99% 
confidence level for half om 12 trial phases. The upper curve shows the limit before any planetary signal is remove to the RV time 
series. The sharp decrease in detection sensitivity around the period 5.3 day is caused by the RV signal of Gl 581b. The lower curve 
shows the limit after the best 4-planet Keplerian fit has been subtracted. The sharp decrease in sensitivity around the period 2 day is 
due to sampling. Venus and Mars criterion delineate the habitable zone, shown in blue. The vertical yellow dashed line marks the 
duration of the survey. 
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Fig. 15. Survey sensitivity derived from the combined phase-averaged detection limits on individual stars. Iso-contours are shown 
for 1, 10, 20, 30, 40, 50, 60, 70, 80 and 90 stars. Planet detected or confirmed by our survey are reported by red circles and labeled 
by their names. 
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Table 4. Test for variability 
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4 


2.0 


1.0 


0.956 


0.9 


0.813 


-0.073 


0.9 


0.998 


0.535 


G1406 


3 


5.7 


5.7 


0.745 


3.9 


0.141 


-0.054 


3.7 


0.995 


0.663 


G1413.1 


17 


1.1 


3.0 


1.4 10""* 


93.1 


< 10 ^ 


0.206 


73.9 


0.008 


0.001 


G1433 


50 


0.8 


3.3 


< 10"^ 


985 


< 10"^ 


0.885 


897 


2.4 10"^ 


0.055 


G1438 


12 


1.0 


3.0 


6.8 10""* 


73.2 


< 10 ^ 


-0.656 


58.7 


0.118 


0.039 


G1447 


6 


0.9 


1.3 


0.370 


12.2 


0.032 


0.889 


10.1 


0.754 


0.415 


G1465 


15 


1.7 


2.2 


0.269 


32.4 


0.004 


-0.811 


26.7 


0.055 


0.153 


G1479 


58 


0.9 


4.1 


<io-^ 


1272 


<io-'^ 


-0.240 


1269 


1.000 


0.137 


LHS337 


8 


3.6 


3.3 


0.721 


9.1 


0.243 


1.159 


7.0 


0.327 


0.188 


G1480.1 


8 


3.3 


1.9 


0.960 


3.6 


0.822 


0.422 


3.2 


0.716 


0.286 


G1486 


4 


2.2 


2.7 


0.580 


4.9 


0.180 


0.740 


3.7 


0.824 


0.229 
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Table 4. continued. 



Name 


N 






P{F) 


A constant 


P(y^ ) 

vt constant^ 


Slope 


slope 




FAP 








[m/sl 








rm/s/vrl 






G1514 


8 


0.6 


1.5 


0.025 


50.6 


1.1 10"^ 


-0.383 


46.6 


0.888 


0.368 


G1526 


29 


0.6 


2.9 


< 10-^ 


887 


< 10-^ 


0.104 


881 


1.000 


0.893 


G1536 


12 


0.7 


2.7 


1.4 10""* 


165 


< 10-^ 


-1.049 


117 


0.025 


0.054 


G1551 


32 


1.3 


2.1 


0.007 


136 


< 10"^ 


-0.234 


125 


0.006 


0.352 


G1555 


7 


1.6 


3.1 


0.127 


30.5 


3.2 10"^ 


-0.105 


27.1 


0.838 


0.313 




5 


1.0 


4.1 


0.007 


121 


< 10"^ 


-1.704 


109 


0.910 


0.408 


G1581 


121 


1.0 


9.8 


< 10"^ 


1.5 10"^"^ 


< 10"^ 


-0.440 


1.5 10"^"^ 


1.000 


0.021 


G1588 


21 


0.6 


1.1 


0.009 


63.1 


2.3 10"^ 


0.093 


61.9 


0.992 


0.041 


G1618A 


19 


1.0 


5.4 


2.6 10"^ 


543 


< 10"^ 


3.681 


78.1 


< 10"^ 


< 10"^ 


G1628 


23 


0.6 


3.6 


< 10"^ 


646 


< 10"^ 


0.616 


618 


0.622 


0.068 


G1643 


6 


2.1 


3.0 


0.347 


15.8 


0.007 


-1.753 


6.7 


0.082 


0.042 


G1667C 


143 


1.0 


4.0 


< 10"^ 


2984 


< 10"^ 


0.916 


2801 


< 10"^ 


0.002 


G1674 


44 


0.6 


6.8 


< 10-^ 


6588 


< 10-^ 


-2.166 


6124 


2.0 10""* 


0.142 


G1678 1 A 


11 


0.7 


3.1 


5.5 10"^ 
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229 
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G1876 


52 


0.9 


120 


< io-» 


1.1 10+^ 


< 10-' 


14.304 


1.1 10+^ 


0.058 


0.018 


G1877 


43 


1.3 


4.0 


< io-» 


882 


< lo-"* 


0.020 


881 


1.000 


0.885 


G1880 


8 


0.7 


2.4 


0.004 


107 


< 10" 


-1.270 


27.9 


0.001 


0.004 


G1887 


75 


0.7 


4.3 


< io-» 


4422 


< 10-' 


1.489 


4061 


<io-^ 


0.057 


LHS543 


7 


2.4 


2.9 


0.466 


11.3 


0.079 


1.105 


7.3 


0.206 


0.094 


G1908 


33 


0.6 


1.8 


2.9 10-» 


385 


< 10-" 


-0.352 
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0.125 


0.444 
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7 
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Fig. 16. Periodograms for RV time series with more than 6 measurements. 
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Fig. 16. Periodograms for RV time series (continued). 
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. 16. Periodograms for RV time series (continued). 
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. 16. Periodograms for RV time series (continued). 
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Fig. 16. Periodograms for RV time series (continued). 
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Fig. 16. Periodograms for RV time series (continued). 
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. 16. Periodograms for RV time series (continued). 
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Fig. 16. Periodograms for RV time series (continued). 
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Fig. 16. Periodograms for RV time series (continued). 
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Fig. 16. Periodograms for RV time series (continued). 
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Name 


Period 
[day] 


K 
[m/s] 


^ sin 


FAP 


FAP2 


Gil 


71.9 


1.50 


330 


0.340 


0.236 


G112 


2.25 


4.69 


0.0 


0.936 


0.991 


LHS1134 


3.77 


17.43 


0.1 


0.935 


0.320 


G154.1 


6.22 


4.73 


6.8 


0.731 


0.634 


G187 


16.2 


1.53 


11.4 


0.339 


0.331 


G1105B 


4.04 


5.20 


14.7 


0.713 


0.674 


CD-44-836A 


5.78 


4.43 


0.2 


0.804 


0.879 


LHS1481 


2.05 


6.74 


0.1 


0.535 


0.447 


LP771-95A 


4.95 


16.30 


1.3 


0.985 


0.717 


LHS1513 


28.8 


5.64 


0.0 


0.989 


0.854 


GJ1057 


5.31 


15.23 


0.7 


0.967 


0.958 


G1145 


5.74 


3.69 


0.0 


0.875 


0.354 


G1176 


8.78 


4.71 


1342 


< 10-3 


0.002 


LHS1723 


2.49 


4.44 


0.0 


0.763 


0.581 


LHS1731 


3.29 


4.31 


0.1 


0.819 


0.398 


G1191 


17.8 


2.45 


196 


0.223 


0.181 


G1203 


3.62 


5.31 


0.5 


0.972 


0.965 


G1205 


32.8 


4.47 


2951 


< 10-3 


< 10-3 


G1213 


4.52 


8.42 


0.3 


0.995 


0.532 


G1229 


2.24 


2.19 


46.9 


0.786 


0.761 


HIP31293 


17.9 


3.19 


0.7 


0.678 


0.487 


HIP31292 


3.46 


4.77 


0.0 


0.440 


0.125 


G1250B 


3.79 


56.06 


9.8 


0.791 


0.894 


G1273 


425 


4.56 


406 


< 10-3 


< 10-3 


LHS1935 


8.27 


3.89 


0.1 


0.915 


0.788 


G1285 


2.24 


160.77 


90.8 


0.292 


0.416 


G1299 


2.34 


7.86 


0.4 


0.635 


0.426 


G1300 


5.28 


6.88 


76.3 


0.346 


0.523 


GJ2066 


2.44 


1.52 


2.0 


0.893 


0.909 


GJ1123 


26.7 


10.12 


0.0 


0.605 


0.159 


G1341 


33.2 


3.26 


90.1 


0.200 


0.319 


GJ1125 


4.57 


217.41 


636 


0.435 


0.192 


G1357 


13.5 


4.77 


0.1 


0.849 


0.863 


G1358 


26.0 


12.74 


403 


< 10-3 


0.002 


G1367 


15.3 


2.77 


40.3 


0.193 


0.899 


G1382 


2.88 


4.07 


843 


0.606 


0.527 


G1388 


2.23 


29.88 


2594 


< 10-3 


< 10-3 


G1393 


37.5 


2.40 


140 


0.099 


0.206 


G1413.1 


11.8 


2.17 


29.5 


0.743 


0.586 


G1433 


7.36 


3.49 


393 


< 10-3 


< 10-3 


G1438 


24.1 


3.65 


18.1 


0.819 


0.916 


G1447 


3.09 


2.01 


0.1 


0.981 


0.861 


G1465 


2.14 


2.81 


7.2 


0.651 


0.753 


G1479 


23.1 


4.30 


637 


< 10-3 


< 10-3 


LHS337 


2.53 


5.67 


0.1 


0.267 


0.908 


G1480.1 


3.41 


3.39 


0.1 


0.801 


0.868 


G1514 


15.2 


2.85 


0.6 


0.353 


0.892 


G1526 


49.5 


3.73 


352 


0.058 


0.059 


G1536 


3.35 


3.81 


33.1 


0.819 


0.992 


G1551 


2.04 


2.53 


65.4 


0.489 


0.202 


G1555 


6.67 


5.42 


0.7 


0.835 


0.825 


G1569A 


9.37 


6.25 


0.3 


0.887 


0.718 


G1581 


5.37 


12.62 


2091 


< 10-3 


< 10-3 


G1588 


4.49 


1.08 


25.8 


0.499 


0.382 


G1618A 


10000 


16.45 


79.4 


0.002 


0.418 


G1628 


67.3 


4.32 


201 


0.033 


0.063 


G1643 


3.55 


10.28 


0.0 


0.938 


0.697 


G1667C 


7.20 


3.96 


1396 


< 10-3 


< 10-3 


G1674 


4.70 


9.27 


1528 


< 10-3 


< 10-3 
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Table 5. continued. 



Name 


Period 
[day] 


K 
[m/s] 


2 


FAP 


FAP2 


G1678.1A 


73.8 


6.66 


11.1 


0, 


.171 


0, 


.239 


G1680 


10000 


27.59 


91.9 


0, 


.001 


0, 


.452 


G1682 


392 


4.99 


8.9 


0, 


.852 


0, 


.978 


G1686 


3.51 


4.19 


0.1 


0, 


.816 


0, 


.624 


G1693 


14.5 


2.35 


0.1 


0, 


.096 


0, 


.899 


G1699 


10000 


63.66 


56.3 


0, 


.280 


0, 


.699 


G1701 


22.2 


4.69 


15.2 


0, 


.205 


0, 


.421 


LTT7434 


7.96 


24.94 


6.7 


0, 


.228 


0, 


.412 


G1729 


2.89 


32.52 


48.1 


0, 


.438 


0, 


.508 


G1752A 


51.3 


7.86 


16.7 


0, 


.032 


0, 


.150 


Lrl/j4 


I.JJ 


D.J / 


U.J 


0, 


.907 


0, 


.399 


LrJ iZJo 


J.OZ 


i i.OJ 


U.U 


0, 


.057 


0, 


.665 


\JJ iZJO 


Z.U/ 


1 1 9Q 


U.U 


0, 


.923 


0, 


.942 


i-/no J Jo J 


Z.UD 






0, 


.944 


0, 


.997 


i^r o iO-DU 


9 Q7 




u.u 


0, 


.232 


0, 


.305 


LrloJz 




1 < /1Q 

I J Ay 


jjZ 


< 


10 ^ 


< 


10 ^ 


G184o 


10.7 


0.32 


O O /C 

336 


< 


10"^ 


< 


10"^ 


G1849 


2000 


30.68 


229 


< 


10 ^ 


< 


10 ^ 


GJ1265 


3.24 


23.22 


0.0 


0, 


.546 


0, 


.287 


Ulo /D 


ol.O 


1 no Q o 


Z.l ID 


< 


10-^ 


< 


10-^ 


G1877 


33.0 


4.61 


477 


0, 


.529 


0, 


.540 


G1880 


10.1 


4.02 


3.7 


0, 


.838 


0, 


.643 


G1887 


13.7 


3.15 


3094 


0, 


.084 


0, 


.020 


LHS543 


38.0 


4.79 


0.1 


0, 


.907 


0, 


.611 


G1908 


9.39 


1.67 


228 


0, 


.937 


0, 


.887 


LTT9759 


2.78 


6.78 


0.8 


0, 


.851 


0, 


.956 
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Fig. 17. Periodograms for RV time series with 1st Keplerian signal removed. 
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Periodograms for RV time series with 1st Keplerian signal removed (continued). 
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Table 6. Test for periodicity after subtraction of the best keplerian fit 



Name 


Period 


K 


2 

^ sin 


FAP 




[day] 


[m/s] 






G1176 


40.1 


4.19 


537 


< 10"^ 


G1205 


39.3 


1.91 


1524 


< 10"^ 


G1273 


17.3 


1.58 


199 


0.038 


G1358 


2.08 


4.21 


158 


0.278 


G1433 


5.76 


1.61 


256 


0.415 


G1479 


11.3 


2.07 


341 


0.040 


G1581 


12.9 


3.64 


946 


< 10-^ 


G1667C 


28.1 


1.72 


1096 


0.012 


G1674 


26.2 


3.69 


540 


< 10-3 


G1832 


35.7 


1.68 


324 


0.016 


G1846 


30.8 


4.73 


142 


0.014 


G1849 


19.6 


1.61 


132 


0.999 


G1876 


30.1 


47.52 


2.6 10-^4 


< 10-3 
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Table 8. Model comparison based on false-alarm probabilities (FAP). Tested models are composed of either a constant, 1,2, 3 or 4 planets. Their 
FAP are in boldface when the more complex model is found to be a statistically significant improvement over the simpler model. 



1x1 o 

iM diiie 






v^sie 


i pidnei 


z pidneis J pidneis 


GU 


1 planet 


8.3 


67.2 








2 planets 


4.1 




94.6 






3 planets 


2.7 


61.2 


9.9 




G154.1 


1 planet 


0.3 


26.1 






Lrlo / 


1 planet 




OJ. / 








1 planet 


U. / 








vJli / U 


i pidllCL 


9S Q 


n 1 








9 nlanpts 


9.2 




0.2 












1 






pidneis 






^ 7 
J. / 


JZ.J DO. J 


niiQi 

VJll !/ 1 


1 t^lllTK^t 

1 piaiiCL 


7 Q 


S9 9 








9 nlanpts 


1.8 




19.0 




G1205 


1 nlanpt 


22.8 


< 10"^ 








9 r^lQTi/=»i"c 
Z pidllCLS 


1 9 ZL 




^ 10-3 






J pidneis 


O.U 


"V lU 


^ in-3 

■V J-U 






4 planets 
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^ 1 A— 3 
< lU 


^ 1 A— 3 


^ 1 A— 3 'TO c 
< lU /J.J 




1 planet 


1. / 










1 planet 


1 n /I 
iU.4 


^ 1 A— 3 
< lU 








2 planets 






Jo.l 






3 planets 


2.6 


< 10 ^ 


46.1 






4 planets 


1.4 


< 10"^ 


34.6 


76.8 87.5 


G1300 


1 planet 


3.0 


21.3 








2 planets 


0.9 




69.6 




Lrl341 


1 planet 


A A 

4.4 


oU.j 






G1358 


1 planet 


17.4 


< 10"^ 








2 planets 


2.9 




2.6 




Lrlio/ 


1 planet 


1 o 
1.0 


30.0 






CjIJoz 


1 planet 


ZU.O 


19.1 








2 planets 


/.J 




ZS'.O 




ul3oo 


1 planet 


OO 1 

OZ.I 


^ 1 A— 3 








2 planets 


/I rv o 

40.8 




1 

67.1 






3 planets 


18.9 


< 10 ^ 


70.3 




G1393 


1 planet 


5.5 


13.4 








2 planets 


1.6 




61.5 




G1413.1 


1 planet 


1.3 


19.4 






G1433 


1 planet 


11.6 


< 10"^ 








2 planets 


6.8 




o o o 

88.2 






3 planets 


3.8 


1.3 


o 

65.8 






4 planets 


1.7 


0.3 


15.3 


/I /I 1 O A O 

AAA 84.3 


G1438 


1 planet 


0.7 


27.6 
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1 planet 


0.6 


86.1 
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1 planet 


12.5 


< 10 ^ 








2 planets 
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13.7 






3 planets 
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1 planet 
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2 planets 
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3 planets 
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< 10-3 






4 planets 


3.4 


< 10-3 


< 10-3 


< 10-3 0.5 


G1588 


1 planet 


0.7 


20.9 






G1618A 


1 planet 


4.3 


< 10-3 






G1628 


1 planet 


4.7 


0.8 






G1667C 


1 planet 


10.0 


< 10-3 
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Table 8. continued. 



NT Jimp 




red 




1 T^ljinpf 
J. jjidiid 




2 planets 


6.3 




< 10-3 




3 planets 


4.7 


< 10-3 


< 10-3 




4 planets 


3.5 


< 10-3 


< 10-3 


G1674 


1 planet 


41.5 


< 10-3 






2 planets 


15.3 




6.0 




3 planets 


7.3 


< 10-3 


1.0 


G1680 


1 planet 


4.7 


2.0 




G1682 


1 planet 


0.4 


89.2 




G1699 


1 planet 


2.6 


84.1 




G1701 


1 planet 


1.1 


10.2 




G1752A 


1 planet 


0.8 


8.6 




G1832 


1 planet 


13.6 


< 10-3 






2 planets 


6.0 




9.8 




3 planets 


2.8 


< 10-3 


1.9 




4 planets 


1.5 


< 10 3 


34.7 


G1846 


1 planet 


14.1 


< 10 3 






2 planets 


4.0 




12.0 


G1849 


1 planet 


7.4 


< 10-3 






2 planets 


2.8 




66.9 


G1876 


1 planet 


2569.0 


< 10-3 






2 planets 


102.3 




< 10-3 




3 planets 


44.5 


< 10-3 


< 10-3 




4 planets 


12.2 


< 10-3 


< 10-3 


G1877 


1 planet 


10.3 


3.0 






2 planets 


5.0 




94.5 




3 planets 


2.0 


0.9 


60.0 


G1887 


1 planet 


36.2 


0.3 






2 planets 


18.9 




2.0 




3 planets 


10.8 


< 10-3 


0.3 




4 planets 


5.4 


< 10-3 


< 10-3 


G1908 


1 planet 


5.3 


33.8 






2 planets 


1.8 




50.4 



2 planets 3 planets 



14.8 



36.6 



52.6 



75.5 



2.6 



0.2 



< 10 



1-3 



8.4 
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Table 9. Model comparison based on false-alarm probabilities (FAP). Tested models are composed of either a constant, 1,2, 3 or 4 planets. Their 
FAP are in boldface when the more complex model is found to be a statistically significant improvement over the simpler model. 



Name 


Model 




Drift 


1 pl.+dr. 


Gil 


1 pl.+dr. 


8.1 


75.3 






2pl.+dr. 


4.3 




92.3 




3 pl.+dr. 


2.0 


37.3 


13.3 


G154.1 


1 pl.+dr. 


0.3 


65.9 




G187 


1 pl.+dr. 


0.9 


71.1 




G1105B 


1 pl.+dr. 


0.7 


57.0 




G1176 


1 pl.+dr. 


26.6 


0.4 






2 pl.+dr. 


9.4 




0.7 




3 pl.+dr. 


5.5 


< 10"^ 


1.1 




4 pl.+dr. 


2.4 


< 10-^ 


< 10-3 


G1191 


1 pl.+dr. 


6.5 


11.8 






2 pl.+dr. 


2.2 




81.1 


G1205 


1 pl.+dr. 


18.7 


< 10-^ 






2 pl.+dr. 


10.1 




< 10-3 




3 pl.+dr. 


6.3 


< 10-^ 


< 10-3 




4 pl.+dr. 


5.3 


< 10-^ 


< 10-3 


G1229 


1 pl.+dr. 


1.5 


91.4 




G1273 


1 pl.+dr. 


10.8 


4.0 






2 pl.+dr. 


4.9 




71.6 




3 pl.+dr. 


2.5 


0.3 


59.0 




4 pl.+dr. 


1.4 


1.1 


68.8 


G1300 


1 pl.+dr. 


2.8 


34.1 






2 pl.+dr. 


0.7 




62.2 


G1341 


1 pl.+dr. 


2.4 


14.9 




G1358 


1 pl.+dr. 


18.3 


< 10-3 






2 pl.+dr. 


4.8 




50.0 


G1367 


1 pl.+dr. 


0.8 


6.8 




G1382 


1 pl.+dr. 


20.7 


17.9 






2 pl.+dr. 


7.1 




56.6 


G1388 


1 pl.+dr. 


85.4 


< 10-3 






2 pl.+dr. 


42.3 




51.6 




3 pl.+dr. 


16.8 


< 10-3 


31.0 


G1393 


1 pl.+dr. 


4.8 


9.6 






2 pl.+dr. 


1.8 




98.1 


G1413.1 


1 pl.+dr. 


1.5 


57.8 




G1433 


1 pl.+dr. 


9.8 


< 10-3 






2 pl.+dr. 


4.9 




44.0 




3 pl.+dr. 


3.7 


< 10-3 


6.2 




4 pl.+dr. 


1.6 


0.1 


32.4 


G1438 


1 pl.+dr. 


0.8 


46.9 




G1465 


1 pl.+dr. 


0.6 


77.0 




G1479 


1 pl.+dr. 


12.8 


< 10-3 






2 pl.+dr. 


7.3 




41.7 




3 pl.+dr. 


4.2 


< 10-3 


14.9 




4 pl.+dr. 


2.5 


0.3 


35.1 


G1526 


1 pl.+dr. 


11.2 


2.6 






2 pl.+dr. 


2.9 




58.0 


G1536 


1 pl.+dr. 


1.8 


89.9 




G1551 


1 pl.+dr. 


1.5 


6.4 






2 pl.+dr. 


0.8 




99.1 


G1581 


1 pl.+dr. 


18.9 


< 10-3 






2 pl.+dr. 


8.6 




< 10-3 




3 pl.+dr. 


5.7 


< 10-3 


< 10-3 




4 pl.+dr. 


4.5 


< 10-3 


< 10-3 


G1588 


1 pl.+dr. 


0.8 


33.7 




G1618A 


1 pl.+dr. 


0.8 


16.7 




G1628 


1 pl.+dr. 


5.0 


0.8 




G1667C 


1 pl.+dr. 


8.8 


< 10-3 





1 pl.+dr. 2 pl.+dr. 3 pl.+dr. 



39.6 27.9 



2.3 59.9 



98.0 74.7 



83.0 



87.3 97.6 



< 10-3 96.6 
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Table 9. continued. 



Name 


Model 




Drift 


1 pl.+dr. 




z pi.+ar. 


J.O 




n 9 




J pl.+clr. 


'I o 




^ 1 A— 3 
< 10 




4 pl.+dr. 


1. 1 
J.i 






G1674 


1 pl.+dr. 


32.5 


< 10 ^ 






z pl.+dr. 


lZ.4 




1 /.z 




3 pl.+dr. 


3.5 


< 10 ^ 


. -i A— 3 
< 10 ^ 


LrlooO 


1 pl.+dr. 


1.0 


43.3 




LjIOoZ 


1 pl.+dr. 


n /L 


8/1 ^ 
o4.0 






1 r\l -Lrlr 

i pi.+ui. 


n Q 


J) J). J 




vjl /Ul 


1 pl.+dr. 


U.o 


/ J. 1 








7 


94 Q 






1 pl.+dr. 


i+.U 








9 nl -\-c\r 


6 9 




lO.U 




J pi.+ur. 


9 8 
z.o 








4 r\] -I-Hr 


1 Q 




7 Q 


G1846 


1 pl.+dr. 


14.6 


0.3 






2 pl.+dr. 


2.8 




1.8 




1 nl -l-rlr 


7 7 


< 10"^ 






2 pl.+dr. 


2.3 




35.1 


vJlO /D 


1 pl.+dr. 


ZH-JJ.Z 








2 pl.+dr. 


y /.I 




^ 1 A— 3 
< 10 




J pl.+dr. 


JZ.J 


< 10 


^ 1 A— 3 
< 10 




4 pl.+dr. 


10. z 


^ 1 A— 3 

< 10 


^ 1 A— 3 
< 10 


Lrio / / 


1 pl.+dr. 




z. / 






2 pl.+dr. 


A A 

4.4 




on c\ 
oV.O 




J pi.+ur. 


i .J 


1 Q 8 
Vy .O 


oo. i 


G1887 


1 pl.+dr. 


34.2 


0.5 






2 pl.+dr. 


17.0 




1.4 




3 pl.+dr. 


9.0 


< 10-3 


< 10-3 




4 pl.+dr. 


4.9 


< 10-3 


< 10-3 


G1908 


1 pl.+dr. 


5.2 


28.8 






2 pl.+dr. 


1.6 




48.4 



2 pl.+dr. 3 pl.+dr. 



0.3 



94.1 



63.4 



97.0 



5.4 



0.1 



5.2 



30.2 
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Fig. 18. Conservative detection limits on m sini for time-series 
with more than 4 measurements. Planets above the limit are ex- 
cluded, with a 99% confidence level, for all 12 trial phases. Some 
panels appear with 2 curves : the upper one is the detection lim- 
its before any model is subtracted and the bottom one is for the 
residuals around a chosen model (composed of planets, linear 
drifts and/or simple sine function). See Sect.[6]for details. 
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Fig. 19. Phase-averaged detection limits on m sini for time- 
series with more than 4 measurements. Planets above the limit 
are statistically excluded, with a 99% confidence level, for half 
the 12 trial phases. Some panels appear with 2 curves : the upper 
one is the detection limits before any model is subtracted and 
the bottom one is for the residuals around a chosen model (com- 
posed of planets, linear drifts and/or simple sine function). See 
Sect. [6] for details. 



Appendix A: Comparison with published time 
series 

A.1. Variability 

Compared to other published time-series, we measured lower 
dispersions for all M dwarfs but Gl 846 and known planet-host 
stars. Gl 1 is not found variable in E06 and Z09 but their dis- 
persion is limited by a higher photon noise (~ 2.6 m/s, against 
(Te = 1.9 m/s in our case). We report variability for G1229 but 
at a level of <1.9 m/s while the variability reported in E06 and 
Z09 implies a jitter of 3.9 - 4.7 m/s. The slightly lower disper- 
sion we observe for G1357 (cTe = 3.2 m/s) against 3.7 and 5.3 
m/s for E06 and Z09, respectively, might not be significant given 
our low number of observations (6). For G1551, we measured a 
dispersion only slightly lower (3.3 against 3.6 m/s). We observe 
significantly lower variability for G1682 (1.8 against 3.6 m/s) 
and Gl 699 (1.7 against 3.4 and 3.3 m/s), and a higher dispersion 
for Gl 846 (5.6 against 3.0 m/s). Although diff'erent time-spans, 
epochs of observations and activity levels at those epochs could 
explain diff'erent dispersions f or ind ividual stars (as it is certainly 
the case for Gl 846 - see Sect.|5.2|), the fact we measure a lower 



slope 


1 nis paper 


Z09 


, 1 • rr . 

(r-aiitcrcncc 




[m/s/yr] 


[m/s/yr] 




Gil 


+0.332 + 0.212 


-0.204 + 0.305 


0.09 


G1229 


-0.257 + 0.296 


+ 1.410 + 0.269 


4.30 


G1357 


-1.682 + 0.710 


+0.273 + 0.305 


2.77 


G1551 


-0.234 + 0.162 


+0.715 + 0.135 


4.44 


G1682 


+0.685 + 0.490 


+2.395 + 0.562 


2.54 


G1699 


-3.043 + 0.646 


-0.697 + 0.133 


3.73 



Table A.l. Linear trends for the time series of stars common to 
2009) Z09) and this paper. The fourth col- 



Zechmeister et al. 



umn reports the significance of the diff'erence, expressed in cr 
(and corresponding to the overlap of 2 gaussian distributions, 
evaluated with Monte Carlo trials). 



dispersion for most comparison stars most likely reflects the bet- 
ter performances of the Harps spectrograph. 

A.2. Trends 

Like this paper, Z09 reports non- significant slopes for G1357 
and G1682 and significant slopes for Gil, Gl 551 and G1699 
(although in our case G1699 is attributed a significant trend by 
the F-test only). Nonetheless, the slopes reported for G1551 
and G1699 seem diff'erent and they moreover found a signif- 
icant trend for G1229 whereas we do not. Time series have 
also been published for Gl 832 and Gl 849 as they were singled 
out from their sample t o report an orbiting planet (Bai ley et al. 
2009 [ Butler et al.|2006i ). For both stars, the planet reflex motion 
clearly dominates the radial velocity signal so we discard them 
from a quantitative comparison. In Table |A.1| we compare the 
slopes of linear fits to the time series in Z09 and to those of this 
paper. We note that most often the significant diff'erences reflect 
a signal more complex than a simple linear drift. 



A. 3. Periodicity 

Among stars with identified periodicity in RV data, G1832, 
Gl 849 and Gl 876 have time series published to report on de- 
tected planets. The periodicities we have found for those three 
stars are similar to their planets' orbital periods. Only Gl 876d is 
not detected with our automated procedure because one has to 
do a full N-body integration to subtract properly the signal in- 
duced by planets 'b' and 'c'. Besides known planet hosts, Z09 
also report on absence of periodicities for Gl 229, Gl 357, Gl 433 
and G1682, and significant periodicities for G1551 and G1699. 
Our results and Z09 are therefore in contradiction for three stars 
: G1433, G1551 and G1699. We noted in §[5l]that, for G1433, 
the RVs reported by Z09 and in this paper are not incompati- 
ble provided that the merged data set is fitted by a model com- 
posed of 1 planet plus a quadratic drift. In the case of G1551, 
the ~l-year periodicity Z09 and jEndl & Kurster] (|2008 ) have at- 
tributed the signal to an alias o f a low frequency signal with the 
typical 1-year sampling. After [Endl & Kurster| ([T008), the low 
frequency signal is believed to be caused by a cluster of points 
that are both blue- shifted and with a higher Ha index compared 
to other points in the time series. This putative activity signal 
might not be seen in our time series because it counts only 24 
measurements, against 229 in Z09. Finally, the periodicity found 
for Gl 699 is also attributed to activity, with a clear counterpart in 
Ha filling factor. Again, if that activity signal is not seen in our 
time series, it is likely because it counts only 22 measurements, 
against 226 for Z09. 



